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ABSTRACT 



Context. Cosmological structure formation offers insights on all of the universe's components: baryonic matter, dark matter, and, 
notably, dark energy. The mass function of galaxy clusters at high redshifts is a particularly useful probe to learn about the history of 
structure formation and constrain cosmological parameters. 

Aims. We aim at deriving reliable masses for a high-redshift, high-luminosity sample of clusters of galaxies selected from the 400d sur- 
vey of X-ray selected clusters. Weak gravitational lensing allows a comparison of the mass estimates derived from X-rays, forming 
an independent test. Here, we will focus on a particular object, CL0030+2618 at z = 0.50. 

Methods. Using deep imaging in three passbands with the Megacam instrument at MMT, we show that Megacam is well-suited for 
measuring gravitational shear, i.e. the shapes of faint galaxies. A catalogue of background galaxies is constructed by analysing pho- 
tometric properties of galaxies in the g'r'i' bands. 

Results. We detect the weak lensing signal of CL0030+2618 at 5.8cr significance, using the aperture mass technique. Furthermore, 
we find significant tangential alignment of galaxies out to ~ 10' or > 2 r 2 oo distance from the cluster centre. The weak lensing centre 
of CL0030+2618 agrees with several X-ray measurements and the position of the brightest cluster galaxy. Finally, we infer a weak 
lensing virial mass of M 200 = 7.2+ 3 «+ 2 - 3 x 10 14 M o for CL0030+2618. 

Conclusions. Despite complications by a tentative foreground galaxy group in the line of sight, the X-ray and weak lensing estimates 
for CL0030+2618 are in remarkable agreement. 



Key words. Galaxies: clusters: general - 
lensing - X-rays: galaxies: clusters 



1. Introduction 



Galaxies: clusters: individuals: CL0030+2618 - Cosmology: observations - Gravitational 



The mass function n(M, z) of galaxy clusters is a sensitive probe 
to both cosmic expansio n and evolution of structure by gravita- 
tional collapse (cf. e.g., Rosati et al. 2002; Voit 2005; Sch ueckerl 
2005b . Therefore, mass functions derived from statistically well- 
understood cluster samples can be and are frequently used to de- 
termine cosmological parameters like D. m , the total matter den- 
sity of the universe in terms of the critical density, and <x 8 , the 
dispersion of the matter density contrast. In addition, measure- 
ments of the mass function at different redshifts have the poten- 
tial to constrain the p ossible evolution of the dark energy compo - 
nent of the universe dPeacock et al.ll2006l : lAlbrecht et al.l F2006). 
expressed by the value and change with time of the equation-of- 
state parameter 

Because the abundance and mass function of clusters are 
sensitive functions of these cosmological parameters, they have 



been studied intensively both theoretically (tress & schechtei 
1974; ISheth&Tormenll 19991; Uenkins et all 120011; iTinker et al 
20081) and observationally. 



* Observations reported here were obtained at the MMT Observatory, 
a joint facility of the Smithsonian Institution and the University of 
Arizona. Our MMT observations were supported in part by a donation 
from the F. H. Levinson Fund of the Peninsula Community Foundation 
to the University of Virginia. In addition, MMT observations used for 
this project were granted by the Smithsonian Astrophysical Observatory 
and by NOAO, through the Telescope System Instrumentation Program 
(TSIP). TSIP is funded by NSF. 



Several methods to measure cluster masses have been em- 
ployed to determine their mass function. Assuming hydrostatic 
equilibrium of the intracluster medium (ICM), the mass of a 
cluster can be computed, once its X-ray gas density and temper- 
ature profiles are known. If the quality of the X-ray data does not 
allow the determination of profiles for individual clusters, for in- 
stance at high redshift, X-ray scaling relations of the X-ray lumi- 
nosity (Lx-M tot ; e.g., Reiprich & Bohringer 2002; Mantz ~et~al1 
2008), temperature (Tx-M tol ), gas mass (M gas -M tot ) or Yx _ 
TxMgas with the total mass are used as proxies dVikhlinin et al.l 
l2009al) . Simultaneously constraining cosmological para meters 
and X-ray cluster scaling relations, iMantz et al. (2009) found 
wde = —1.01 ±0.20 for the dark energy equation-of-state param- 
eter, compiling data from a large sample of galaxy clusters. 

Weak gravitational lensing provides a completely indepen- 
dent probe of a cluster's mass, as it is sensitive to baryonic and 
dark matter alike, not relying on assumptions of the thermody- 
namic state of the gas. The fact that sources of systematic errors 
in the lensing and X-ray methods are unrelated opens the pos- 
sibility to compare and cross-calibrate X-ray and weak lensing 
masses. 
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Several studies have been recently undert aken, compar- 
ing X-ray and weak lensing cluster observables: iDahld (2006) 
found the weak lensing mass to scale with X-ray luminosity 
like oc /^04*0-4 6 an( j cons t ram a combination of O m and erg. 
lHoekstral (2007) established a proportionality between the weak 
lensing mass M w i within the radius r25oo inside which the den- 
sity exceeds the critical density by a factor 250 and with 
an exp onent a — 1.34^ ^. For the same radius, Ma hdavi et al.l 
(20083 quoted a ratio Mx/M w] = 1.03 + 0.07 along with a de- 
crease towards smaller r adii. Directly aiming at the ratio of weak 
lensing to X-ray mass, Zh ang et al.l (|2008) measured at a ra- 
dius of M w i/Mx = 1.09 + 0.08, for a radius rspp, and later con - 
firmed this value and t he trend with radius (Zhang e t al.ll20 10h . 
ICorless & Kingl d2009h investigated the role of the mass estima- 
tor on the systematics of the weak lensing mass function. 

In order to make progress, it is particularly important to 
determine the masses of more clusters to a high accuracy, es- 
pecially at high (z > 0.3) redshifts. To date, only a few studies 
have been undertaken in this regime which provides the strongest 
leverage on structure formation and is thus crucial for tackling 
the problem of dark energy. 

In the redshift range 0.3 < z < 0.8, for two reasons, weak 
gravitational lensing becomes more advantageous over meth- 
ods relying on the X-ray emission from the ICM with increas- 
ing redshift. First, X-ray temperature profiles are becoming in- 
creasingly difficult to determine with cluster redshift. Second, 
the fraction of clusters undergoing a merger increases with z, 
in accordan ce with the paradigm of hierarchical structure for- 
mation (e.g. ICohn & Whitdl20 051 rendering the assumption of 
hydrostatic equilibrium more problematic at higher z. 

With this paper, we report the first results of the largest 
weak lensing follow-up of an X-ray selected cluster sample at 
high redshifts. 



2. Observations 

2.1. The 4006 Survey and Cosmological Sample 

The 400 Square Degree Galaxy Cluster Survey (abbreviated as 
400d) comprises all clusters of galaxies detected serendipitously 
in an analysis of (n early) all suitable ROSAT PSPC pointings 
dBurenin et alj|2007l) . 

The survey's name derives from the total area of 397° 2 on 
the sky covered by these pointings. The sample is flux-limited, 
using a threshold of 1.4 x 10~ 13 erg s _1 cm -2 in the 0.5-2.0 keV 
band. 

The analysis o f the ROSAT data is described in detail in 
Bure nin et alJ d2007l) . All clusters in the 400d sample are con- 
firmed by the identification of galaxy overdensities in optical 
images. Their redshifts have been determined by optical spec- 
troscopy of sample galaxies. 

The final 400d catalogue contains 242 objects in the red- 
shift range 0.0032 < z < 0.888. In order to be able to accu- 
rately constrain the mass function of galaxy clusters at z ~ 0.5, 
specifications of the 400d sample were devised such that the 
cluster catalogue provides a sample of "typical" clusters in the 
0.3 < z < 0.8 range. 

This work is based on the 400d Cosmological Sample, 
a carefully selected subsample of high-redshift and X-ray lu- 
minous 400d c l usters. It has been defined and published in 
Vikh linin et ail (12009a. Table 1) and comprises 36 clusters. 
This cosmological or high-redshift sample is drawn from the 
400d catalogue by selecting all clusters both at a redshift z > 



0.350, as given bv iBurenin eTaT] (120071) . and with a ROSAT lu- 
minosity exceeding 

£x,min= 4.8 xl0 43 (l+z) 18 ergs" 1 . (1) 

IBurenin et al.l d2007l) assume a ACDM cosmology with Hubble 
parameter h = 0.72 and density parameters O m = 0.30 and £2a = 
0.70 for matter and dark energy, respectively. For consistency, 
the same cosmology is adopted throughout this paper. 

The X-ray luminosity limit in Eq. (Q} has been chosen in 
such a way that it, given an Lx-^x-scaling relation valid for 
low-z clusters, corresponds to a mass of « 10 14 Mq. This renders 
our cosmological sample nearly mass-limited. 

All clusters in the cosmological sample have been reob- 
served with Chandra within the framework of the Chandra 
Cluster Cosmology Project (CCCP) and constraints on cosmo- 
logical pa rameters derived from these X-ray data have been pub- 
lished in I Vikhlinin et all (l2009allbl) . 

With this paper, we report first results of a weak lensing 
follow-up survey of the galaxy clusters in the 400d cosmological 
sample. Here, we focus on one particular object, CL0030+2618, 
which, as we will see below, represents an exceptionally inter- 
esting case. Also, we demonstrate in detail the methods we use 
for data reduction and analysis as this is the first weak lensing 
study using Megacam at MMT. 

The clus ter CL0030+2618 is listed at a redshift of z = 
0.500 both in IBurenin et all d2007l) (designation BVH 002) as 
well as in the precurs or to the 400d, termed the 160d survey 
dVikhUninetalJll998l) as VMF 001 It w as first identified as 
a cluster of galaxies by Boy le et al.l d 1997b conducting a spec- 
troscopic follow-up to Rosat observations in the visual wave- 
length range. These authors assigned the desig nation CRSS 
J0030.5+2618 and measured a redshift of z = 0.5 16. Brandt et al. 
(2000) observed the field of CL0030+2618 with Chandra dur- 
ing its calibration phase, studying fain t hard X-ray sources in 
the vicinity of the cluster. iHorner et al.l d2008l) confirm the red- 
shift of z = 0.500 for the cluster with their designation WARP 
J0030.5+2618 in their X-ray selected survey of ROSAT clus- 
ters, but point out a possible contamination of the X-ray signal 
by a line-of-sight structure at the lower redshift of z~ 021 . 

Addition al Chandra observation s were conducted as part 
of the CCCP dVikhlinin et alJl2009allbh . Its X-ray emission as 
detected by Rosat is centred at aj2ooo = 00 h 30 m 33?6, feooo = 
+26°1 8'16". The analysis of Chandra data bv IVikhlinin et al.l 
(2009a) results in a luminosity in the 0.5 . . . 2.0 keV-band of 
L x = 1.57 x 10 44 ergs -1 and a ICM temperature of k B T x = 
(5.63 + 1.13)ke V. CL0030+2618has b een classified as a possi- 
bly merging bv lVikhlinin eTaT] (|2009a) based on its X-ray mor- 
phology. 

As CL0030+2618 has not been studied with deep optical 
imaging before, and there are no observations with large op- 
tical telescopes in the major public archives, we present the 
first such study of this cluster. (The SEGUE observations used 
in Sect. 13.2.21 for cross-calibration have some overlap with our 
Megacam imaging south of CL0030+2618, but do not contain 
the object itself.) 

2.2. The Megacam Instrument at MMT 

The observations discussed in this work have been obtained us- 
ing the Megacam 36-chip camera at the 6.5 m MMT telescope, 
located at Fred Lawrence Whipple O bservatory on Mt. Hopkins, 
Arizona dMcLeod et al.ll2000l 120061) . Megacam is a wide-field 
imaging instrument with a field-of-view of ~ 24' x 24', resulting 
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from a mosaic of 4x9 CCDs, each consisting of 2048x4608 pix- 
els which corresponds to a very small pixel size of 0.08" px" 1 . 
Each chip has two read-out circuits and amplifiers, each reading 
out half a chip (cf. Sect. IA.ll in the Appendix). The gaps between 
the chips measure 6" in the direction corresponding to declina- 
tion using the default derotation and 33", 5", and 33" in the 
direction associated with right ascension. We use Megacam in 
the default 2x2 binning mode. 

A system of u'g'r'i'z' filters, similar to but subtly dif- 
ferent from their na mesakes in the Sloan Digital Sky Survey 
dFukugita et al l ll996h is used for Megacam. The relations be- 
tween the Megacam and SDSS filter systems are described in 
detail in Sect. lA.5l and visualised in Fig. IA.3I 

None of the previous studies with Megacam (e.g.: 
iHartman et al.l2008l:IWalsh et al.120 08) is related to gravitational 
lensing. Thus, in this paper we will show that Megacam indeed 
is suitable for weak lensing work. 

2.3. Observing Strategy 

In principle, the small distortions of background sources which 
we want to measure are achromatic. In practice, however, the op- 
timal passband for weak lensing observations is determined by 
the signal-to-noise ratio which can be obtained in a given amount 
of time and depends on seeing and instrumental throughput. To 
maximise the number of high signal-to-noise, background galax- 
ies whose shapes can be determined reliably for a given exposure 
time, we choose the r'-band as the default lensing band. Aiming 
at a limiting magnitude of r' lim « 26 for r exp as 3 h, we obtain a suf- 
ficient number of high-quality shape sources (n ga i > 15 arcmuT 2 ) 
in the final catalogue^ 

Lensing effects depend on the relative distanc es between 
source and deflector dBartelmann & Schneiderl200lL Chap. 4.3). 
Ideally, we would like to determine a photometric red shift es 
timate for each galaxy in our lensing catalogue (e.g.: Bemte: 



2^jBolzonell^taf]p002 IWolfet alj200 H lllbert et afl2W 
Hildebrandt et al.l2008l) . However, this is observationally expen- 
sive as deep imaging in > 5 passbands is necessary to obtain 
accurate photometric redshifts. 

On the other hand, using only one filter (the lensing band) 
and a simple magnitude cut for a rough separation of background 
from foreground galaxies needs a minimum of observing time 
but neglects the galaxies' intrinsic distribution in magnitude. We 
are following an intermediary approach here, using three filters 
from which we construct colour-colour-diagrammes of the de- 
tected galaxies and use this information to achieve a more ac- 
curate background selection than using the simplistic magni- 
tude cut. This method has success fully been applied to weak - 
lensing galaxy clus t er data by e.g . Clowe & Schneider! J2002); 
lBrada£ et al] (120051) : iKausch et all (120071) . Megacam' s g' and if 
passbands straddle the Balmer break, the most distinctive feature 
in an elliptical galaxy's optical spectrum in the redshift range 
Z ~ 0.5 in which we are interested. Therefore, we use the g'r'i' 
filters and resulting colours to identify foreground and cluster 
objects in our catalogues. 

In order to obtain a high level of homogeneity in data qual- 
ity over the field-of-view despite the gaps between Megacam's 
chips, we stack dithered exposures. Our dither pattern consists 
of 5 x 5 positions in a square array with 40" distance between 
neighbouring points, inclined by 10° with respect to the right as- 
cension axis on which the chips normally are aligned. We find 



1 Eq. 10 gives our definition of "limiting magnitude", while the ac- 
tual values for CL0030+2618 are listed in Tabled 



this pattern to be robust against missing frames (exposures which 
couldn't be used in the final stack for whatever reason). 

2.4. The Data 

The data presented in this paper have been collected in five 
nights distributed over two observing runs on October 6th 
and 7th, 2004 and October 30th, October 31st, and November 
1st of 2005. In these two observing runs, a total of four 
400d Cosmological Sample clusters have been observed. In the 
first phase of data reduction, the so-called run processing, these 
data have been processed in a consistent fashion. A weak lens- 
ing analysis of the three clusters besides CL0030+2618, namely 
CL0159+0030, CL0230+1836, and CL0809+2811, will be the 
topic of an upcoming paper in this series. In this work, we 
make use of part of these data for the photometric calibration 
of CL0030+2618 images (see Sect. EXT). 

In Table [JJ we give the exposure times, seeings, and re- 
lated information for the final image stacks on which we base 
all further analysis. Fig. [JJ shows a three-colour composite im- 
age prepared from the stacked Megacam g'r'i' observations of 
CL0030+2618. 



3. Outline of Data Reduction 

The reduction of the data presented in this paper relies on the 
THELI pipeline originally designed and tested on observations 
obtained using th e Wide-Field Ima ger (WFI) on ESO's La Silla 
2.2 m telescope (Erben etalJ 120051) . While the reduction fol- 
lows the procedure detailed in lErben et al.l d2005l) . some impor- 
tant changes had to be made to adapt the THELI pipeline to 
work on MMT Megacam data. In the following discussion of the 
subsequent individual reduction routines, special emphasis has 
been given to those adaptations arising from the fact that MMT 
Megacam is a "new" camera with a small field-of-view per chip 
(325"X 164" instead of 853" x 379" for MegaPrime at CFHT, or 
a factor 1 /6 in field-of-view), using a larger telescope. 

The THELI pipeline distinguishes two stages of data reduc- 
tion: run processing and set processing. During run processing, 
the first phase, all frames taken during an observation run in a 
particular filter are treated in the same way. Run processing com- 
prises the removal of instrumental signatures, e.g. de-biasing and 
flatfielding. In set processing the data are re-ordered according 
to their celestial coordinates rather than their date of observa- 
tion. Astrometric and photometric calibration lead to a resulting 
"coadded" (stacked) image for each set. In Fig. [2] we present the 
coadded r'-band image of CL0030+2618, overlaid with its final 
masks in the left and its weight image in the right part of the plot. 

3.1. Coaddition "Post Production" 

The final stage of the data reduction is to mask problematic re- 
gions in the coadded images, applying the methods presented in 
Dietr ich et al.l (2007). By subdividing the image into grid cells of 
a suitable size and counting SExtractor detections within those, 
we identify regions whose source density strongly deviate from 
the average as well as those with large gradients in source den- 
sity. This method not only detects the image's borders but also 
masks, effectively, zones of increased background close to bright 
stars, galaxies or defects. 

In a similar fashion, we mask bright and possibly satu- 
rated stars which are likely to introduce spurious objects in cat- 
alogues created with SExtractor. We place a mask at each posi- 




Fig. 1. Three-colour composite of CL0030+2618, prepared from the Megacam g'r'i' coadded images. The main image shows a cut- 
out of the central region of CL0030+2618, with an edge length of »255" (1600px), corresponding to the virial radius of w 1.5 Mpc 
at the cluster redshift of z- 0.50. North is up and East is to the left. The tentative luminous arcs near the galaxies Gl and G3 (Table|2] 
and Sect. 15. 4t are emphasized in the two smaller, zoomed (40" x 40") images. 

Table 1. Specifications of the coadded images for CL0030+2618 



Filter 


Observation Dates 


CM 


Seeing 


Calib. Method 


mim 


r' 


2004-10-06/7 


6600 


0.82 


stellar colours 


25.9 


g' 


2005-10-30/1, 2005-11-01 


7950 


0.87 


SDSS standards 


26.8 


i' 


2005-10-31 


5700 


1.03 


SDSS standards 


25.1 



tion of such sources as drawn from the USNO Bl catalogue. The 
method in which the size of the mask is scaled acc ording to the 
star's magnitude has been described in some detail in lErben et alj 
(2009). A small number of objects per field which are missing in 
the USNO Bl catalogue has to be masked manually while masks 
around catalogue positions where no source can be found have 

to be removedQ 

Conforming with Hildebrand t et alJ d2005l) . we compute the 
limiting magnitudes in the coadded images for a 5cr-detection in 
a 2" aperture as 

OTi im = Zf - 2.5 log (5 V-^pixO-sky) (2) 

with Zf being the photometric zeropoint in the respective filter 
/, /Vpi x the number of pixels within the 2" aperture and cr^y the 
RMS sky -background variation measured from the image. 



2 This last manual step can be largely avoided by also automatically 
masking objects drawn from the Hubble Space Teles cope Guide Star 
Catalog (GSC), as demonstrated bv lErben et"aTJ d2009h . 



Obtaining accurate colours for objects from CCD images 
is not as trivial as it might seem. Next to the photometric cali- 
bration (see Sect. 13.21 1. aperture effects have to be taken into ac- 
count. Our approach is to measure SExtractor isophotal (ISO) 
magnitudes from seeing-equalised images in our three bands. 
We perform a simplistic PSF matching based on the assumption 
of Ga ussian PSFs as it has been described in Hildebrandt et al. 
d2007h : The width of the filter with which to convolve the k-th 
image is given as 

°" filter.k = y]crl orst - °~l (3) 

with o\ and cr worst being the widths of the best fitting Gaussians 
to the PSFs measured from the A:-th and the worst seeing image. 

3.2. Photometric Calibration 
3.2.1. Calibration Pipeline 

The photometric calibra tion of our data is large ly based on the 
method developed by Hildebran dt et alJ d2006l) but using AB 
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Fig. 2. Left: The coadded r'-band image of CL0030+2618, overlaid with its final masks. The target cluster is located in the frame 
centre. Small square masks cover regions masked because of their source counts strongly deviating from the average in the field 
(Sect. 13.11 ). The elongated masks enclose tracks of slowly moving objects (asteroids) which had to be identified on the image by 
visual inspection. The small octagonal masks are saturated stars found from the USNO Bl catalogue and manually. Right: The r'- 
band weight image of CL0030+2618. Pixels lying in the interior of the chips get significantly higher weights than those which fall 
into an intra-chip gap in some of the dithered exposures. 



magnitudes and SPSS-like filters together w ith the SDSS Data 
Release Six (lAdelman-McCarthv et alJl2008l DR6) as our cali- 
bration catalogue. 

As mentioned in Sect. |24] the CL0030+2618 data we 
present here have been observed together with three other clus- 
ters from the 400d sample. Two of these, CLO 159+0030 and 
CL0809+2811, are situated within the SP SS PR6 footprint. In 
addition to these, numerous IStetsonl {2000) standard fields were 
observed along with the science data of which some provide fur- 
ther SPSS photometry. 

The MMT/Megacam filter system is based on the SPSS one 
but not identical to it (see Fig. IA.3b : the comparison between the 
systems is detailed in Sect. IA.5l in the Appendix. Therefore, re- 
lations between instrumental magnitudes and calibrated magni- 
tudes in the SPSS system have to take into account colour terms. 

In order to determine the photometric solution, we use 
SExtractor to draw catalogues from all s cience and stan- 
dard f rames having SPSS overlap. Using the Hildeb randt et al.l 
(2006) pipeline, we then match these catalogues with a photo- 
metric catalogue assembled from the SPSS archives, serving as 
indirect photometric standards. 

The relation between Megacam instrumental magnitudes 
rai nst and catalogue magnitudes wsdss for a filter / can be fit- 
ted as a linear function of airmass a and a first-order expansion 
with respect to the colour index, simultaneously: 

Wf.inst - w fjSDSS =y3f(m f>S DSS - «f,SDSs) + Ha + Z f (4) 

where mf SDSS - ?«f',SDSS is a general colour index w.r.t. another 
filter /', /3f the corresponding colour term, and Zf the photomet- 



ric zeropoint in which we are mainly interested. For the fit, we 
select objects of intermediate magnitude which are neither satu- 
rated nor show a too lar ge scatter in mj nst given a certain twsdss- 
Following the model of Hilde brandt et al.1 d2006l) . we account for 
the variable photometric quality of our data by fitting /3f, jf, and 
Zf simultaneously under the best conditions, fixing jf in interme- 
diate, and fixing jf and /3f in even poorer conditions. The fixed 
extinctions and colour terms are set to default values which are 
discussed in Sect. IA.5l in the Appendix. 

3.2.2. Colour-Colour-Diagrammes 

Comparing the zeropoints for different nights and fields, we 
conclude that the nights on which the r'-band observations of 
CL0030+2618 were performed, were not entirely photomet- 
ric but show a thin, uniform cirrus. Therefore, an indirect re- 
calibration method is needed here. To this end, we fitted the 
position in the r 1 — i' vs. g' — r 1 colour-colour-diagramme of the 
stars identified in the CL0030+2618 field to those found in two 
other, fully calibrated, galaxy cluster fields, CL0159+0030 and 
CL0809+2811. 

In the left panel of Fig. [3] we plot the g' - r' versus r' - V 
colours of stars identified in these two fields and compa re them 
to theo retical spectra of main-sequence stars from the iPicklesI 
( 1998) spectral library to find good agreement between both the 
two observed sequences and the predicted stellar colours. 

As we find reliable absolute photometric calibrations for 
the g'- and /'-bands of CL0030+2618, the location of the stellar 
main sequence for this field is determined up to a shift along the 
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• CL0030+2618fjeld 


• CL0030+2618 MMT 






: * CL01 59+0030 field 
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■ CL0030+2618SEGUE 
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♦ Pickles (1998} 




2 


7 ♦ Pickles (1998) * t ^ : 




2 


1 




■J. 


1 





♦ 

♦♦ 

4 


♦ 
♦ 
♦ 

* 
♦ 

♦♦ 

4 





1 






-1 
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Fig. 3. Photometric calibration by stellar colours: Left panel: plotted here are the g'-r' vs. r'-i' colours of sources identified as stars 
in three galaxy cluster fields observed with Megacam. For two of these fields, CL0159+0030 (upward triangles) and CL0809+281 1 
(downward triangles), absolute photometric calibration with SDSS standards could be performed. For CL0030+2618, results for 
recalibrated r'-band zeropoints ar e shown (dots; details see main text) The colours in all three fields agree with the colours 
of main sequence stars from the iPicklesi dl998l) spectral library (diamonds). Right panel: The g' - r' vs. r' - i' colours of 
stars in the Megacam images of CL0030+2618 (dots) w hich could also be identified in the partially overlapping SEGUE strip 
dNewberg & SloanD igital Skv Su rvey Collaborationll2003l) and shown here as squares are both consistent with each other as well 
as with the lPicklesi Q.1998) colours (diamonds). Each pair of measurements of one individual source is connected. 



main diagonal of the g'-r' versus r'-i' diagramme, corresponding 
to the r' zeropoint. We fix this parameter by shifting the main 
sequence of CL0030+2 618 on top of the other observed main 
sequences as well as the lPicklesi (119981) sequence. We go in steps 
of 0.05 magnitudes, assuming this to be the best accuracy we 
can reach adopting this rather qualitative method and settle for 
the best-fitting test value (see Table[T]i. The dots in Fig.[3]show 
the best match with the CL0159+0030 and CL0809+28 1 1 stellar 
colours obtained by the re-calibration of the CL0030+2618 r'- 
band. 



After the photometric calibration, we became 
aware of a field observed in the SEGUE projec t 
dNewberg & Sloan Digital Skv Survey Collaboration] 120031) 
using the SDSS telescope and filter system which became 
publicly available along with the Sixth Data Release of SDSS 
dAdelman-McCarthv et alJ [2008) and partially overlaps with 
the CL0030+2618 Megacam observations. Thus, we are able 
to directly validate the indirect calibration by comparing the 
colours of stars in the overlapping region. The right panel of 
Fig. [3] shows the good agreement bet ween th e two i ndependent 
photometric measurements and the Pickles (1998) templates 



from which we conclude that our calibration holds to a good 
quality. 

For comparison we also calibrated the CL0030+2618 r'- 
band by comparing its source counts to the ones in the 
CL0159+0030 and CL0809+2811 fields for the same filter, but 
discard this calibration as we find a discrepancy of the r esulting 
main s equence in g'-r' versus r'-i' with the theoretical IPicklesi 
(1998) models mentioned earlier. 

4. The Shear Signal 

Gravitational lensing leads to a distortion of images of distant 
sources by tidal gravitational fields of intervening masses. Here, 
we describe the method to measure this shear while referring to 
ISchneiderl d2006l) for the basic concepts and notation. 

4.1. KSB analysis 

The analysis of the weak lensing data is based on the 
iKaiser et alJ dl995l KSB) algorithm. The reduction pipeline we 
use was adapted from the "TS" implementati on presented in 
iHevmans et al.l (I2006I) and explored further in Schrabbac k et alJ 
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0[px] 

Fig. 4. The distribution of sources in the apparent size - mag- 
nitude - space. Plotted are the SExtractor magnitudes r' AmQ 
against the half-light radii ft of all sources in the "reliable" cat- 
alogue. The stellar locus is visible prominently and is defined 
for further analysis by 2.55 px < ft < 2.95 px and 16.75 mag < 
r AUTO < ma 8 ( trie § re y sna ded area). The sources to the 
right and lower sides of the thick lines r' Al]T0 < 16.75 mag and 
ft > 2.8 px for sources fainter than 22.5 mag we treat as un- 
saturated galaxies. In addition, the magnitude limits mbright an d 
Wfaint defined in Sect. IB. 21 are given as thin lines for illustrative 
purposes. 



d2007l) : lHartlap et alJ (120091) . Its basic concepts are already out- 
lined in lErben et al.l d2001l) . Thus, in this section, we will focus 
more on the properties of our data than on the methods them- 
selves since they are well documented in the above references. 

The KSB algorithm confronts the problem of reconstruct- 
ing the shear signal from measured galactic ellipticities; there- 
fore, it has to disentangle the shear y from the intrinsic ellip- 
ticitie^l of the galaxies and from PSF effects. The simultaneous 
effects of shear dilution by the PSF and the superposition of the 
incoming ellipticity with the anisotropic PSF component can be 
isolated by tracing the shapes of sources we can identify as stars, 
bearing neither intrinsic ellipticity nor lensing shear. 

The complete correction yielding a direct, and ideally unbi- 
ased, estimator s of the (reduced) shear g exerted on a galaxy in 
our catalogue reads: 

^ = (p g Uh-p|;((p sm *);i4)] (5) 

Here, small Greek indices denote either of the two components 
of the complex ellipticity. Quantities with asterisks are measured 
from stellar sources. The ep are the ellipticities of galaxies as 
measured from the input "shape" image while the e* g are the el- 
lipticities of stars tracing the PSF (plotted in the left two panels 
of Fig. |5]). The 2x2 matrices P s and P sm give the transformation 
of ellipticities under the influences of gravitational shear fields 
and the presence of an (anisotropic) PSF, respectively. See, e.g., 
ISchneiderl d2006l) on how these quantities are determined from 
higher-order moments of the measured brightness distributions. 



3 In this study, we adopt the following definition of ellipticity: if r< 1 
is an ellipse's axis ratio its ellipticity is described by a two-component 
(polar) quantity e with \e\ = (1 - r 2 )/(l + r 2 ) which we represent as a 
complex number with Cartesian components e = e\ + ie2- 



4.2. The KSB and Galaxy Shape Catalogues 

Catalogues are created from the images using the SExtractor 
double detection mode: Sources are identified on the lensing 
band image in its original seeing. Photometric quantities (fluxes, 
magnitudes) are determined at these coordinates from the mea- 
surement images in the three bands g'r'i' convolved to the worst 
seeing (found in the /'-band). 

The photometric properties determined from the three 
bands are merged into one catalogue based on the detection im- 
age. From those catalogues problematic sources are removed. 
Those are sources near the boundaries of the field-of-view or 
blended with other sources, as well as objects whose flux radii 
do not fall in the range ft* <ft g < 10 px, for which KSB works 
safely, with ft* the angular size of unsaturated stars. For the re- 
maining objects, the shapes can now be determined. 

Note that the KSB catalogue presented in Fig.|4]and all cat- 
alogues discussed hereafter only contain those objects for which 
a half- light radius ft could be determined by our implementation 
of the Er ben et all (1200 ll) method. Objects for which the mea- 
surements on the (noisy) data yield negative fluxes, semi major 
axes, or second-order brightness moments or which lie to close 
to the image border are removed from the catalogue, reducing its 
size by ~3 %. 

Figure [4] shows the distribution of the sources in the "reli- 
able" catalogue in apparent size - magnitude space. The promi- 
nent stellar locus enables us to define a sample of stars by 
ft* . <ft< ft* and r'* < r'. lrrr , < r'* with ft* . = 2.55 px, 

mm max mm AUTO max mm " ' 

^max-2.95 px, /*• = 16.75 mag, and r^ ax =22.5 mag (the shaded 
area in Fig. HJ from which the PSF anisotropy et in Eq. (0 is 
determined. 

In the transition to the galaxy shape catalogue, we regard as 
unsaturated galaxies all objects ^ UT0 > (i.e. fainter than the 
brightest unsaturated point sources) and more extended than ft> 

#max for 'auto < r max or 2-8 px for ^uto > r max' respectively. 
The latter is justified by the fact that while for bright sources 
it is easy to distinguish galaxies from point sources, there is a 
significant population of faint galaxies for which a very small 
radius is measured by the SExtractor algorithm. Thus, we relax 
the radius criterion by 5 % for sources fainter than r^ ax . 

However, we notice that among those small objects there is 
a population of faint stars, not distinguishable from poorly re- 
solved galaxies using an apparent size - magnitude diagramme 
alone and resulting in a dilution of the lensing signal compared 
to a perfect star - galaxy distinction. Our decision to neverthe- 
less include these small sources into our catalogue is based on 
the resulting higher cluster signal as compared to a more con- 
servative criterion (e.g. ft/ft* mx < 1-10 for the galaxies fainter 
than r^ ax ). We call "galaxy shape catalogue" the list of objects 
that both pass this galaxy selection and the cuts for signal qual- 
ity discussed in Sect. 14.61 This important catalogue yields the 
final "lensing catalogue" by means of the background selection 
discussed in Sect. 14.41 

4.3. PSF Anisotropy of Megacam 

In the KSB pipeline, we fit a model e cot '*(x,y) in the pixel co- 
ordinates x and y to the measured ellipticities e* of stars such 
that the residual anisotropies e anl '* = e* — e cor '* of stellar im- 
ages should effectively be zero. Figure|5]shows the effect of PSF 
anisotropy correction: The raw ellipticities of the tracing stars 
as presented in the left two panels are modelled by a polynomial 
g cor _ Yj"Zq pux k y e defined globally over the whole field-of- 
view. The best-fitting solution for the case n = 5 we adopt here is 
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Fig. 5. Correction of PSF anisotropy of the CL0030+2618 r' band as used in the analysis. The upper panel shows the distribution of 
the ellipticity components of the stars identified in the field, together with the numeric values of their dispersions <t\ 2 : = a{e\ ,2). 
The "whisker plots" in the lower panel show the size and orientation of PSF anisotropy as they vary as a function of the spatial 
coordinates x and y. On the left, the situation before correction, i.e. the ellipticities as measured in the stars are depicted. The middle 
two plots give the fit by a global fifth order polynomial in x and y. Residuals of this correction are presented in the plots on the right. 



shown in the middle panels of Fig. [5] while the residual elliptici- 
ties of the stars e' dm are displayed in the panels to the right. 

Simultaneously aiming at reducing both the mean (el m '*) 
of the residual ellipticities and their dispersion cr(e™ *) we find 
that a polynomial order as high as n — 5 is necessary to effec- 
tively correct for the distinctive quadrupolar pattern in the spa- 
tial distribution of the "raw" stellar ellipticities (see lower left 
and middle panels of Fig. [5]). There is no obvious relation be- 
tween the zones of preferred orientation of the PSF ellipticity 
in Fig. |5] and the 4x9 chip detector layout of Megacam. See 
Sect. l4.6l for further details. 

For stacking in the lensing band, we select only those 
frames which show moderate PSF ellipticity in the first place 
(see Sect. I A. 4l in the Appendix for details). Thus, we ensure the 
images used for lensing analysis to be isotropic to a high de- 
gree even before any corrections are applied. By stacking im- 
ages in which the PSF anisotropy is different in magnitude and 
orientation (cf. Figs. IA.ll and lA.2K we further reduce the ellip- 
ticity owing to the imaging system. The total amount of PSF 
anisotropy present in our Megacam data is small: Before cor- 
rection, we measure <<?i) = 1.77 x 1(T 5 , (e 2 ) = -4.03 x 1(T 3 , 
<|e|> = 1.10 x 10~ 2 , and cr(|e|) = 6.19 x 10~ 3 , reducing after 
the correction to (e'f) = -5.60 x 10~ 7 , (ef) = -2.60 x 10~ 5 , 
<|e ani |> = 4.34 x 10~ 3 , and cr(|e ani |) = 2.90 x 10~ 3 . Note that the 
very small averages for the individual components result from 



partial cancellation of anisotropies from different parts of the 
field-of-view. Thus, MMT/Megacam shows a similar degree of 
PSF anisotropy as other instruments from which lensing sig- 
nals have be en measured successfu lly, e.g. MegaPrime/Megacam 
on CFHT dSemboloni et aD l2006h or Subaru's SuprimeCam 
(Oka be & Umetsul 120081) . The latter authors measured, as an 
RMS average of seven galaxy cluster fields, (ei) = 1.41 X 10~ 2 , 
<e 2 ) = 1 -63 X 10~ 2 , and <x(|e|) = 2.32 x 10~ 2 before correction with 
larger values for the anisotropy components but a simpler spatial 
pattern. 



Although we find small-scale changes in the PSF ellipticity 
which have to be modelled by a polynomial of relatively high 
order, the more important point is that the PSF anisotropy varies 
smoothly as a function of the position on the detector surface in 
every individual exposure, showing a simpler pattern than Fig. [5] 
(See Fig. lA.2l for example exposures at both small and large val- 
ues of overall PSF anisotropy induced by the tracking behaviour 
of MMT) Consequently, it can be modelled by a smooth func- 
tion which is a necessary prerequisite for using the instrument 
with the current weak lensing analysis pipelines. Thus, we have 
shown that weak lensing work is feasible using MMT Megacam. 
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Fig. 6. The effect of the polynomial correction for the PSF 
anisotropy on the ellipticities of galaxies averaged in equally 
populated bins. As a function of the amount of correction e cor 
applied to the components 5-1 (left panels) and 5-2 (right 
panels), we show the raw ellipticities before correction in the 
upper panels and the PSF-corrected ellipticities s in the lower 
panels. The bars in the abscissa and ordinate denote the range of 
the bin and the standard deviation of the ellipticity in this bin, 
respectively. 



4.4. Selection of lensed background galaxies 

Before we proceed with the details of our lensing analysis, we 
explain how we arrive from the galaxy shape catalogue at the 
"lensing catalogue" of objects we classify as background galax- 
ies w.r.t. to CL0030+2618. This background selection, as we will 
call it from now on= is based on their g'r'i' photometry. While 
unlensed objects remaining in the catalogue dilute the shear sig- 
nal, rejection of actual background galaxies reduce it as well. 
Note that a sensible foreground removal is especially important 
for relatively distant objects like the 400d Cosmology Sample 
clusters. 

We introduce two free parameters in our analysis: the mag- 
nitude limit mf a i nt below which all fainter galaxies are included in 
the shear catalogue, regardless of their g'-r' and r'-i' colour in- 
dices, and the magnitude m^ght above which all brighter galax- 
ies will be considered foreground objects and discarded. Only 
in the intermediary interval mbright < r' < mf a ; nt the selection of 
galaxies based on their position in the colour-colour-diagramme 
will take place. In these terms, a simple magnitude cut would 
correspond to mbright = Mfaint- We vary these parameters in order 
to optimise the detection of CL0030+2618 and find mbright = 20.0 
and mf a i n t = 22.5. For details of the colour-colour-diagramme 
method, see Sect. IB. 21 The photometric cuts reduce the cata- 
logue size by 6.0 %, leaving us with a lensing catalogue of 
A^cat = 14813 objects, corresponding to a galaxy surface density 
of « C at = 21.2 arcmin~ 2 . 



4.5. Aperture mass and lensing detection 

The weak lensing analysis we conduct is a two-step process. 
First, we confirm the presence of a cluster signal by constructing 
aperture mass maps of the field which will provide us with a po- 
sition for the cluster centre and the corresponding significance. 
In the second step, building on this position for CL0030+2618, 



the tangential shear profile can be determined and fitted, leading 
to the determination of the cluster mass. 

More precisely, we use the so-called S -statistics, corre- 
sponding to the signal-to-noise ratio of the aperture mass esti- 
mator which for any given centre C is a weighted sum over the 
tangential ellipticities of all lensing catalogue galaxies within a 
circular apertu re of radius ou t . The estimator can be written an- 
alytically as (ISchneiderl ll996): 



(6) 



where s t j denotes the measured shear component tangential with 
respect to the centre for the galaxy at position As filter func- 
tion Q(x := Ifl j —flcl /flout), we apply the hyperbolic tangent filter 
introduced by Schi rmer et al.l d2007l) : 



<2tanh(*) : 



1 



tanh (x/x c ) 



1 + e a - hx + e 



c+dx 



with the width of the filter determined by x c — 0. 15 and the shape 
of its exponential cut-offs for small and large x given by the de- 
fault values {a, b, c, d}-{6, 150, 47, 50). The S -statistics includes 
as a noise term the intrinsic source ellipticity, calculated from the 
data galaxies as o- e = {e\+s^ 1 I 2 with typical values <x e a:0.38. 

The value of out in Eq. © is also fixed such that it max- 
imises Se oul (6 c ) which strongly depends on the filtering size used. 
Exploring the parameter space spanned by 6 0Ut and the photo- 
metric parameters mbright and mf a ; nt , we find, independent of the 
latter two, the highest S -values with 14' < 8 out < 15'. The be- 
haviour of S as a function of out (at a fixed C ) i s in g ood gen- 
eral agreement with the results of ISchirmer et al.l d2007l) for the 
same filter function QtanhM- Thus, we fix 6 0Ut = 14'5 for the 
further analysis, noting this number's agreement with the size of 
our Megacam images (cf. Fig . [8j . We also tested the influence of 
the parameter x c in the Qtanh filter and find that, with all other 
parameters kept fixed, in the 0.15 < x c < 0.6 interval the max- 
imum S -value changes by less than 0.5 % but decreases more 
steeply for smaller values of x c . 

Applying these parameters and measuring S on a reference 
grid of 60" mesh size, we detect CL0030+2618 at the level of 
5.8cr in a grid cell centred at a distance of 34" from the ROSAT 
position at aj2ooo = 00 h 30 m 33 s .6, <5 J2 ooo = 26°18'16", i.e. smaller 
than the grid resolution. We will investigate further into the clus- 
ter position in Sect. [ 



4.6. Verification of the shear signal 

In this subsection, we summarise the consistency tests per- 
formed on the data to validate the galaxy shape measurements 
giving rise to the shear signal discussed below. 

- Correction of PSF anisotropy: We assess the performance 
of the correction polynomial by analysing the PSF-corrected 
ellipticities e™j 5 of galaxies as a function of the amount 
of correction e™ r that has been applied to them by fitting 
a polynomial to the anisotropy distribution of star images 
(see Sect. I4.3l l. Theoretically, the expected positive corre- 
lation between the uncorrected ellipticities and the correct- 
ing polynomial should be removed and e™](e cor ) thus scatter 
around zero. We note that most anisotropy is found in the 
6 = 2 component from the beginning (Fig. [6]). This is removed 
in the corrected ellipticities, with <e a "j 2 > = -0.0010+0.0010 
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Table 2. Notable galaxies in the field of CL0030+2618. 



Galaxy 


oikooo 


<5j2000 


,./MMT 

'auto 


„/MM'l 

6 AUTO 


,-/MMT 

'auto 


z 


Note 


See Figure 


Gl 


OO h 30 m 34 s .O 


26° 18' 09" 


19.20 


21.14 


18.31 


0.5 16 1 


dominant in CL0030+2618 


Fig. 1 


G2 


00 h 30 m 37.9 


26° 18' 18" 


18.82 


20.23 


18.27 


n.a. 


dominant in foreground group 


Fig. 1 


G3 


OO h 30 m 3e , :3 


26° 19' 20" 


19.46 


20.76 


18.95 


n.a. 


strong lensing feature 


Fig. 1 


G4 


oo^o'^s 


26°20'56" 


17.23 


17.98 


16.94 


0.493 1 


QSO 


Fig. 9 



1 Redshift taken from lBovle et al.l dl997h . 
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Fig. 7. The S -statistics (solid line) as a function of the maximum 
value of the ellipticity estimator s resulting we include in the 
galaxy shape catalogue. The dashed and dash-dotted lines show 
the sizes of the resulting catalogue before and after background 
selection (see Sect. IB.2l i. respectively. 



marginally consistent with zero in the standard deviation. In 
the (5=1 component, we measure a residual anisotropy of 
<e™ j) = -0.0026 + 0.0010 which is one order of magnitude 
smaller than the lensing signal we are about to measure. 
Alternatively to the n — 5 polynomial correction to the entire 
image, we consider a piecewise solution based on the pattern 
of preferred orientation in Fig. [5] Dividing the field into four 
regions at y = 6 1 00 px and at x — 4300 px for y < 6 1 00 px and 
x = 5800px for y > 6100px with a polynomial degree up to 
n — 5 we do not find a significant improvement in terms of 
(e am '*), 0"(e am '*), or e™](e cor ) over the simpler model defined 
over the whole field. 

Maximum shear: Due to the inversion of the noisy ma- 
trix P 8 in Eq. (|5), resulting values for the estimator |e| are 
not bound from above while ellipticities are confined to 
< e < 1 . Thus, attempting to measure weak lensing using 
the KSB method, we need to define an upper limit max(|e|) 
of the shear estimates we consider reliable. We evaluate 
the influence of the choice for max(|e|) on the S -statistics 
(Eq.(6j> by varying it while keeping the other parameters, like 
min (tr P g ), the minimum min (v) of the signal-to-noise ratio 
v of the individual galaxy detection determined by the KSB 
code, and the photometric parameters mbdght and mf a j nt de- 
fined in Sect. IB.2l fixed. In the range 0.6 < max (|e|) < 0.8, we 
find an increasing shear signal due to the higher number of 
galaxies in the catalogues using less restrictive cuts (Fig. [7]). 
For max(|e|)>0.8, we see a sharp decline of the lensing sig- 



nal which we explain as an effect of galaxies entering the 
catalogue whose ellipticity estimate is dominated by noise. 
We fix max(|e|) = 0.8, min(trPS) = 0.1, and min(v) = 4.5 
simultaneously to their given values. We note that, while op- 
timising the S -statistics, this might introduce a bias in the 
mass estimate as a cut in max(|e|) = 0.8 directly affects the 
averaging process yielding the shear. 

Shear calibration: We can account for this bias by scaling 
the shear estimates with a shear calibration factor f such 
that e— »/oe to balance biases like the effect of max (|s|). The 
question how gravitational shear can be measured unbiased 
and precisely has been identified as the cruc ial challenge to 
future weak lensing experiments (see e.g. , iHevmans et al.l 
2006; Mas sev et al.l2 007; Bri dle etalj2009l) . The "TS" KSB 
method employed here has been studied extensively and is 
well understood in many aspects. In order to correct for 
the biased shear measurements, fo und by testing the KSB 
pipeline with the simulated data in Hevma ns et al.l ([2006), 
the shear calib ration factor has been introduced and studie d 
subsequently dSchrabback et alj 120071: lHartlap et ail [2009). 
As pointed out by these authors, the calibration bias depends 
on both the strength of the shear signal under inspection, as 
well as on the details of the implementation and galaxy se- 
lection for the shear catalogue. In the absence of detailed 
shape measurement simulations under c luster lensing condi- 
tions, we chose a fiducial /o = 1 .08 from lHartlap et al.l d2009l) 
and assign an error of <x fo = 0.05 to it, covering a significant 
part of the discussed interval. 

Complementary catalogue: We check the efficacy of the 
set of parameters we adopted by reversing the selection of 
galaxies and calculating the S -statistics from those galax- 
ies excluded in our normal procedure. Reversing the back- 
ground selection, i.e. only keeping those galaxies regarded 
as cluster or foreground sources, we find from 10 s bootstrap 
realisations of the complentary catalogue an aperture mass 
significance of S = -0.83 ± 1.06. From the consistency with 
zero, we conclude these cuts to effectively select the signal- 
carrying galaxies. As the background selection removes / p h = 
^Vcompiem/^Vcat = 6.0 % of the sources in the catalogue, we only 
expect a small bias as / p hScompiem/S'cat ~ -0.8 % resulting 
from the background selection. 



5. The Multi-Wavelength View of CL0030+2618 

5.1. Identifying the BCG of CL0030+2618 

Figure[T]shows two candidates for the brightest cluster galaxy of 
CL0030+2618, galaxies with extended cD-like haloes and sim- 
ilar /'-magnitudes (Tabled. The galaxy Gl, closer to the Rosat 
and C handra cent r es of C L0030+2618, was attributed to a clus- 
ter bv lBovle et al.l d!997l) . measuring a spectroscopic redshift of 
Zgi =0.516, while three of their six spectro-zs are z~ 0.25. 
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Table 3. Colours of prominent galaxies observed in the 
CL0030+2618 field compared to colours computed from 
CWW80 elliptical templates at z = 0.50 and z = 0.25. 



Galaxy 




z 


g'-i' 


g'-r' 


r'-i' 


Gl 




0.516 


2.83 


1.94 


0.89 


G2 




n.a. 


1.96 


1.41 


0.55 


CWW80 Ellz= 


= 0.50 


0.50 


2.782 


1.876 


0.906 


CWW80 Ellz= 


= 0.25 


0.25 


2.098 


1.554 


0.544 



We note that Gl and G2 show different colours in Fig. Q] 
each being similar to their fainter immediate neighbours. As very 
extended sources, Gl and G2 are flagged early-on in the pipeline 
but are included in the raw SExtractor catalogues. Aware of their 
larger uncertainties, we use these magnitudes^ for Gl, G2, and 
two other interesting extended galaxies (Table|2]i. 

The observed g'-r', r'-i', and g'- /'colours are compared 
to the on es predicted for a typ ical BCG at z = 0.50 and z = 0.25, 
using the lColeman et all d 19801 CWW80) elliptical galaxy tem- 
plate (Tablef3]l. Nicely consistent with its spectroscopic redshift, 
we find the colours of Gl to be similar to the z = 0.50 tem- 
plate, while G2's bluer colours resemble the CWW template at 
z = 0.25. We conclude that Gl, located close to the X-ray cen- 
tres, is a member of CL0030+2618, and indeed its BCG. On the 
other hand, G2 can be considered the brightest member of afore- 
ground group at z ~ 0.25 . The existence of such foreground struc- 
ture is corroborated by the broad g'-i' distribution (Fig. IB. II ). Its 
implications are discussed in Sect. l5.2l and l6.3.l1 

5.2. Comparing Centres of CL0030+261 8 

The S -statistics lensing centre We determine the centre of the 
CL0030+2618 lensing signal and its accuracy by bootstrap re- 
sampling of the galaxy catalogue of N cat galaxies used in the 
measurements of the S -statistics. From the basic catalogue we 
draw 10 5 realisations each containing A^ cat sources. For each 
realisation, we determine the S -statistics in the central region 
of 3' x 3' side length (~ IMpc or roughly the virial radius of 
CL0030+2618) using a gridsize of 15" and record the highest 
S -value found on the grid and the grid cell in which it occurs. 

Re-running 10 5 bootstrap realisations of the lensing cata- 
logue with the centre fixed to the lensing centre, we calculate a 
detection significance of S(0 WL ) = 5.94 ± 0.98. 

Weak Lensing Mass Reconstruction In order to get an impres- 
sion of the (total) mass distribution in CL0030+2618, we per - 
form a finite field mass reconstruction ( Seitz & Schneider I2001h . 
This method directly aims at the two-dimensional mass distribu- 
tion k(0) and breaks the mass-sheet degeneracy, i.e. the fact that 
the reduced shear, our observable, is invariant under a transfor- 
mation k — > /c( 1 - A) + A with an arbitrary scalar A, by assuming 
k = along the border of the field. 

The resulting mass map, derived by smoothing the shear 
field with a scale of 2' is shown in Fig. [8] and a zoomed version 
displaying the central region of CL0030+2618 as Fig. [9] The 



4 Here, we use SExtractor AUTO instead of ISO magnitudes, known 
to be more robust at the expense of less accurate colour measurements. 
Nevertheless, we find only small differences between the two apertures, 
allowing for cautious direct comparison. 



thick contours give the surface mass densityjf] Beside the clear 
main peak of CL0030+2618, we find a number of smaller addi- 
tional peaks whose significance we are going to discuss in the 
following section. 

Chandra and X MM -Newton We compare these lensing re- 
sults to detections by two X-ray observatories, Chandra and 
XMM-Newton. For Chandra, we use a surface brightness map 
produce d from the 58 ks ACIS exposure by Vikhlinin et al. 
(2009a) (me dium- thick, blue lines in Figs. [8]and|9j. Using the 
IZhangetalJ feolO) method, we find the flux-weighted Chandra 
centre at 00 h 30 m 34 s .9, 26°18'05", slightly off the flux peak at 
00 h 30 m 34 s .0, 26° 18' 13". 

For XMM-Newton, we show detections in the EPIC-MOS2 
chip, binned in 64 x 64 pixels and smoothed with an adaptive 6<x 
Gaussian kernel. Therefore, the respective contours (medium- 
thick, magenta lines in Figs. [8] and [9]! appear more jagged. 

Lensing and X-ray Centres As can be seen from the cross 
in Fig. [9] the cluster centre determined with the aperture mass 
technique falls within the most significant (k > 0.05) conver- 
gence contour and is, within its lcr error ellipse of 24" x 21", 
in good agreement with the flux-weighted Chandra centre of 
CL0030+2618, separated by 17". Similar, it is consistent with 
the Rosat centre just outside the confidence ellipse and the 
XMM-Newton contours ( |5). All of these cluster centres are, 
in turn, within <20" distance from Gl, the BCG. 

Optical Galaxy Light In addition, we determine the distribution 
of r'-band light from galaxies by adding the fluxes of all un- 
flagged sources in the SExtractor catalogue whose magnitudes 
and flux radii are consistent with the criteria defined for the 
galaxy catalogue in Sect. 14.21 and Fig. [40 We do so for each 
pixel of an auxiliary grid, then smoothing it with a Gaussian of 
60" full-width half-maximum. In Figs.[8]and|9] the r'-band flux 
density is given in isophotal flux units per Megacam pixel, with 
a flux of one corresponding to a r 1 = 25.8 galaxy assigned to 
that pixel (thin red contours). There is, amongst others, a dis- 
cernible r'-band flux peak centred between the galaxies Gl and 
G2 (Fig.|9>. 

5.3. Secondary Peaks 

The shear peak clearly associated with CL0030+2618 is the 
most dominant signal in the Megacam field-of-view, in the lens- 
ing /e-map as well as in the X-rays, which can be seen from the 
XMM-Newton count distribution. In the smoothed r'-band light 
distribution, CL0030+2618 shows up as a significant but not the 
most prominent peak. We have to stress that the background se- 
lection using the m^ght an d m faint parameters optimise the lens- 
ing signal for CL0030+2618, with the likely effect that cluster 
signals at other redshifts and hence with different photometric 
properties will be suppressed. Keeping this in mind, we compare 
secondary peaks in the K-map to apparent galaxy overdensities, 



5 The small surface mass densities, in contrast to the fact that 
CL0030+2618 likely has strong lensing arcs (see Sect. 15.4b . hinting at 
k> 1 locally, is due to smoothing. 

6 In absence of a usable half-light radius & for the more extended 
galaxies, we have to substitute flux radii r g here. Using the observed 
relation between & and r g in our dataset, we consider as galaxies objects 
with r,,>3.5px at 16.75<r^ UTO <22.75 and r„>3.2px at r^ UTO >22.75. 
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Fig. 8. The r'-band image of CL0030+2618, overlaid with r'-band galaxy light (thin, red), Chandra (medium-thick, blue; within the 
smaller square footprint), and XMM-Newton (medium-thin, magenta), and lensing surface mass density contours (thick, green). We 
show X-ray surface brightness levels in multiples of 5 x 10~ 9 cts cm^s^'arcmkT 2 in the 0.5 ... 2. keV band. The r'-band flux density 
contours (thin red lines) start from 0.015 flux units per pixel, in intervals of 0.005 flux units. Lensing convergence contour levels 
were obtained smoothing the shear field y(ff) with a Gaussian filter of 2' width and are linearly spaced in intervals of Ak = 0.01, 
starting at /c = 0.01. XMM-Newton contours show MOS2 counts smoothed by an adaptive Gaussian kernel in logarithmic spacing. 
The labels "PI" to "P10" designate the peaks discussed in Sect. 15.31 



as indicated by the smoothed distribution of r'-band light, and to 
the X-ray detections. 

The galaxy listed as G4 in Table [2] a strong X-ray 
emitter detected with a high signal by both Chandra and 
XMM-Ne wton, is id e ntified as a QSO at redshift z - 
0.493 by iBovle et al. ( 19911) and confirmed to be at z = 
0.492 by ICappi et al] (2001) who found a significant over- 
density of 0.5 ... 2 keV Chandra sources in the vicinity of 
CL0030+2618. Regarding its redshift, it is thus a likely mem- 
ber of CL0030+26 1 8 . 

The Chandra analysis finds two additional sources of ex- 
tended X-ray emission at low surface brightness One of them, 



"PI" in Fig. [8] (see Table [4] for coordinates of this and all 
following peaks) is also detected by XMM-Newton and has 
been identified as a probable high-redshift galaxy cluster by 
iBoschinl d2002l) (his candidate #1 at a J20 oo = 00 h 3 l m 01 s .3, 6 m oo = 
26°26'39") in a deep survey for galaxy clusters using pointed 
Chandra observations. In the k map, contours near the north- 
eastern corner of Megacam's field-of-view extend close to the 
position of this cluster, but their significance near this corner and 
close to the bright star BD+25 65 is doubtful. The Megacam im- 
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Fig. 9. Zoomed version of Fig. [8] showing only the central region of CL0030+2618.The cross gives the position and lcr uncertainty 
of the centre position of the S -peak obtained by bootstrapping (cf. Sect. 15.2b while the star symbols with five and six points denote 
the X-ray centres found with ROSAT and Chandra, respectively. The source on the northern edge on which strong X-ray emission 
is centred is named G4 in Table [2] 



ages show a small grouping of r' « 21 galaxies with similar 
colour in the three-colour composite at the position of "PI". 

The other Chandra peak, "P2" is located near a prominent 
peak in the r'-band light, but with a strong contribution from a 
single bright galaxy within its 60" smoothing radius. It does not 
correspond to a tabulated source in either NEE0 or SIMBAE0 
We do not notice a significant surface mass density from lens- 
ing at this position, but have to stress again that a possible signal 
might have been downweighted by the catalogue selection pro- 
cess. 



7 Not visible in Fig.[8]due to its binning. 

8 NASA-IPAC Extragalatic Database: 
http : / /nedwww . ipac . caltech . edu/ 

9 http : // simbad . u-strasbg . fr/ simbad/ 



Most peaks in the k map, apart from the one associated 
with CL0030+2618, are located at a distance smaller than the 2' 
smoothing scale from the edges of the field, likely due to noise 
amplification by missing information. Amongst them, only the 
second strongest k peak, "P3" seems possibly associated with an 
overdensity of galaxies, but the coverage is insufficient to draw 
further conclusions. 

For a shear peak "P4" close to several Chandra and XMM- 
Newton peaks, there also is an enhancement in r'-band flux, 
while galaxies do not appear concentrated. Likewise, the high 
flux density close to a possible shear peak "P5" seems to be 
caused by a single, bright galaxy. 

On the other hand, we notice agglomerations of galax- 
ies ("P6" to "P8") with a cluster-like or group-like appear- 
ance that show neither X-ray nor lensing signal. For "P7", the 
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Table 4. The additional shear, X-ray, and optical flux peaks dis- 
cussed in Sect.P 



.r Ccuv 


CJ2000 


c 

CJ2000 


UCLCCLCU uy 


P1 


nn h ^ i m n? s 

UU J 1 \)L 


Z.U Z.U JU 


yv-iay, optical 


P2 


00 h 31 m 19 s 


26°25'0" 


X-ray, optical 


P3 


00 h 29 m 31 s 


26° 26' 


shear 


P4 


00 h 31 ra 17 s 


26° 18' 


shear 


P5 


00 h 29 m 49 s 


26 D 15'20" 


shear 


P6 


00 h 30 m 54 s 


26° 23' 


optical 


P7 


00 h 31 ra 10 s 


26°21'15" 


optical 


P8 


00 h 31 m 09 s 


26°13'20" 


optical 


P9 


00 h 31 ra 14 s 


26° 10' 30" 


optical 


P10 


00 h 30 m 51 s 


26° 07' 30" 


optical 



nearby XMM-New ton signal is the distant quasar named 13 by 
iBrandt et al.l yOOO). The two strong r'-band flux overdensities 
"P9" and "P10" in the south-east corner of the Megacam image 
appear to be poor, nearby groups of galaxies. 



5.4. Arc-like Features in CL0030+2618 

We note that, being a massive cluster of galaxies, 
CL0030+2618 is a probable strong gravitational lens, leading 
to the formation of giant arcs. Indeed, we identify two tentative 
strong lensing features in our deep Megacam exposures. The 
first is a very prominent, highly elongated arc ~ 20" west from 
the BCG (Fig. Q). Its centre is at a J20 oo = 00 h 30 m 32 s .7 and 
feooo = 26°18'05"; its length is > 20". The giant arc is not 
circular but apparently bent around a nearby galaxy. 

The second feature possibly due to strong lensing is lo- 
cated near galaxy G3 which appears to be an elliptical. With 
the centre of the tentative arc at crnooo = 00 h 30 m 3G i .5 and 5j2ooo = 
26°19'14", it is bent around the centre of the galaxy forming the 
segment of a circle with ~ 6" radius. Thus, an alternative ex- 
planation might be that the arc-like feature corresponds to a spi- 
ral arm of the close-by galaxy. However, this seems less likely 
given its appearance in the Megacam images. If this arc is due to 
gravitational lensing it is likely to be strongly influenced by the 
gravitational field of the aforementioned galaxy as it is opening 
to the opposite side of the cluster centre. 

Whether these two candidate arcs are indeed strong lensing 
features in CL0030+2618 will have to be confirmed by spec- 
troscopy. 



6. Mass Determination and Discussion 

We analyse the tangential shear profile g t (8) = (s t (8)), i.e., the 
averaged tangential component of e with respect to the weak 
lensing centre of CL0030+2618 found in Sect. I5.2l as a function 
of the separation to this centre. At this point, we also introduce 
the shear calibration factor, fo = 1.08, an empirical correction to 
the shear recovery by our KSB method and catalogue selection 
(cf. Sect. 14.61 ), and the contamination correction factor f\ (ff) we 
will specify in Sect. 16.21 thus replacing s by /o/i (ff)s. First, the 
iNavarro et al.l ([l997, NFW) shear profile will be introduced. 



6.1. The NFW model 

To derive an estimator for the mass of CL0030+2618 from the 
weak lensing data, we fit the tangential shear profile g t (9) with a 



r(z=0.5) [Mpc] 
2 3 4 



fitted f(r)oc2 NFW (r)+const. 




10 



15 



en 



Fig. 10. The fraction of "red sequence-like" galaxies 2.2 < g'—i' < 
3.0 as a function of clustercentric distance before (open symbols) 
and after (filled symbols) background selection. The solid line 
denotes the best fitting sum / cg of a NFW surface mass profile 
and a constant to the latter. We use f\ = f CR + 1 as a correction 
factor for cluster contamination in Sect. 16.31 



NFW profile (e.g. Bar telmannlll996t IWright & Brainerdl feOOO). 
The NFW density profile has two free parameter^! the radius 
r2oo inside which the mean density of matter exceeds the critical 
mass density p c by a factor of 200 and the concentration param- 
eter cnfw from which the characteristic overdensity 6 C can be 
computed. 

The overdensity radius r2oo being an estimator for the clus- 
ter's virial radius, we define as the mass of the cluster the mass 
enclosed within r2oo, given by: 



An 3 
M 2 oo = 200 —p c r 200 . 



The reduced shear observable is: 



<?nfw(«) - 



7nfw(«) 
1 - atnfw(m) 



(8) 



(9) 



where the dimensionless radial distance u = CNFW^d^oo con " 
tains the angular separation 9 and the angular diameter distance 
£>d betwee n lens and observer. The Tn fw(«) and knfw(m) profiles 
are given in Wright Ife Brainerdl (|2000). The critical surface mass 
density 



AnGD A \ D s 



(10) 



depends on and the mean ratio (Dd s /D s ) of angular diameter 
distances between source and observer and source and lens. 



6.2. Contamination by Cluster Galaxies 

In addition to the background selection based on g'—r' and r'-i' 
colours we estimate the remaining fraction of cluster galaxies 
in the catalogue using the g' - i' index. We will use this to de- 
vise a correction factor accounting for the shear dilution by (un- 
sheared) cluster members. As discussed in Sect. IB. ll the colour- 
magnitude diagramme of the CL0030+26 18 field (Fig. IB. Q does 



10 While INavarro et alj i 19971) originally designed their profile as a 
single-parameter model, we follow the usual approach in weak lens- 
ing studies of expressing the NFW profile in terms of two independent 
parameters. 
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Table 5. Properties of the fiducial model combining the param- 
eter values and assumptions going into the NFW modelling. 



Parameter 


Value 


see Sect. 


max(|£|) 


0.8 


1461 


min(v) 


4.3 


14.01 


min(trP g ) 


0.1 




"ibright 


20.0 


m 


'Wfaint 


22.5 


u 


centre 


from S -statistics 


EH 


radial fit range 


0'<e<15' 


16.3.11 


fo 
MO) 

(D is /D s ) 


1.08 ±0.05 
f ww i0) + l 

0.33 ± 0.03 


14.61 

m 

lB.3.21 



not show a clear-cut cluster red sequence, but a broad distribu- 
tion in g' — i', indicating two redshift components. We therefore 
define a wide region 2.2 < g' - i' < 3.0 of possible red sequence 
sources, including galaxies with colours similar to the z = 0.50 
CWW elliptical template but redder than the z = 0.25 one (cf. 
Table [3}. As this definition of "red sequence-like" galaxies is 
meant to encompass all early-type cluster members, it will also 
contain background systems, giving an upper limit for the actual 
contamination in the catalogue. 

Figure [10] shows the fraction of sources 2.2 < g' — i' < 3.0 
in the galaxy catalogue before (open symbols) and after (filled 
symbols) the final cut based on m^ght and mf a i nt has been applied 
as a function of distance to the centre of CL0030+2618 as de- 
termined by lensing (Sect. 15.21 ). Error bars give the propagated 
Poissonian uncertainties in the counts. We note a strong increase 
of the number of "red sequence-like" systems compared to the 
overall number of galaxies towards the cluster centre, indicating 
that a large fraction of those are indeed cluster members. Most 
intriguingly, the background selection seems to remove only few 
of these tentative cluster members, with the fractions before and 
after selection consistent within their mutual uncertainties at all 
radii. This finding can be explained to a large extent by galaxies 
too faint to be removed by the background selection criterion: 
If background selection is extended to the faintest magnitudes 
(?Wfaint = 29), no significant overdensity of "red sequence-like" 
galaxies at the position of CL0030+2618 is detected. Although 
using a different selection method, thi s modes t effec t of back- 
ground selection is in agreement with lHoekstral (120071) . 

By repeating this analysis centred on several random posi- 
tion in our field and not finding a significant increase of the "red 
sequence-like" fraction towards these positions we show that the 
peak around the position of CL0030+2618 is indeed caused by 
concentration of these galaxies towards the cluster centre. 

We find the residual contamination to be well represented 
by the sum f cg (6) = f^ g w (0) + f? g of a NFW surface mass profile 
and a constant ( solid line in Fig. ITOt. We follow the approach of 
Hoekstra (2007) and define a radially dependent factor f\{6) = 
/ C ^ FW (0)+ 1 correcting for the residual contamination. Here we 

take into account only the NFW component / C g FW (6>) of the fit, 
as the offset / c o represents a population of field galaxies, and not 
diluting cluster members. This correction factor scales up the 
shear estimates close to the cluster centre, counterweighing the 
dilution by the larger number of cluster members there. 



6.3. Mass Modelling of CL0030+261 8 
6.3.1 . Fits to the Ellipticity Profile 

In Fig. QTJ there is a discernible positive tangential alignment 
signal extending out to ~10' or ~3.5 Mpc) from the cluster cen- 
tre. (The solid line and dots in all panels give the shear averaged 
in bins of 90" width.) In order to validate that this tangential 
alignment is indeed caused by gravitational shear of a cluster- 
like halo, we fit the NFW reduced shear profile given in Eq. (0 
to the measured shear estimates, probing the range 0' <9< 15'. 

We define a fiducial model using the preferred parameter 
values presented in Table[5] The table also lists references to the 
sections where these values are justified. In order to determine 
r2oo and cnfw, we fit an NFW model to the shear estimates of the 
lensing catalogue galaxies, defined by the parameters above the 
vertical line in Table [5] Parameters below the line do not affect 
the catalogue but influence the relation between shear and cluster 
mass. 

The fitting is done by minimising x 1 using a n IDL imple- 
mentation of the Levenberg-Marquardt algorithm dMorel [l978: 
iMarkward I l2009h and returning r®^ = 1.64 ±0.1 6 Mpc and 
c nfw =2.1±1.1 for the free parameters of the model. Comparing 
the best-fitting NFW model (dashed curve in the upper and mid- 
dle panels of Fig. [TTJ to the data, we find the shear profile to 
be reasonably well-modelled by an NFW profile: we measure 
X 2 l v &o{ — 13404/ 13636 ~ 0.98, assuming an error 

<Tfit=/i(0)a-gt , ( r gt =/ (7- s /V2«0.29 (11) 

for the individual shear estimate. This overall agreement with 
NFW is consistent with shea r profiles of cluster s with compara- 
ble redshift and data quality (Clowe et al. 2006). We discuss the 
NFW parameter values obtained by the fit and the radial range 
over which the NFW fit is valid (the middle and lower panels of 

Fig.OJJiinSect.S] 

Gravitational lensing by a single axially symmetric de- 
flector causes tangential alignment of the resulting ellipticities. 
Thus, the ellipticity cross-component g x = (e x (6)) corresponding 
to a pure curl field around the centre should be consistent with 
zero at all 6. The dotted line and diamonds in the upper panel 
of Figure QTJ show that g x is indeed consistent or nearly consis- 
tent with zero in its error bars in all bins but the innermost 90". 
This feature is, like the general shapes of both g t and g x , quite 
robust against the choice of binning. A tentative explanation for 
the higher g x in the central bin might be additional lensing by 
the foreground mass concentration associated with the z « 0.25 
galaxies (cf. Sect.|5TJ, centred to the East of CL0030+2618. 

To further investigate this hypothesis, we split up the ellip- 
ticity catalogue into an eastern (aj2ooo > ffCLOiBo) and western 
(ffj2000 < Q"cloo3o) subset (with 50.0 % of the galaxies in each) 
and repeat the profile fitting for both of them separately, as the 
influence of a possible perturber at the position of G2 should be 
small compared to the eastern sub-catalogue. In accordance with 
the mass distribution displayed in Fig. |8] in which a higher and 
more extended surface mass density can be found west of the 
centre of CL0030+2618 than east of it, the g t signal is more sig- 
nificant in the sources lying to the West of the cluster than to the 
East. We find ^o,w = L82 ± 0-22 Mpc, c* !WW =2.1 ± 1.2, and 
r 200E = !- 47 ± 0.25 Mpc, c^ FWE = 1.5 ± 1.4. The cross compo- 
nents in the central bins of both subsets are similarly high than 
in the complete catalogue with the eastern half also showing a 
high g x in the second bin. As the values for r2oo from the two 
sub-catalogues are consistent given their uncertainties, we find 
no clear indications for a significant impact of the foreground 
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structure. The inconspicuous lensing signal is consistent with the 
inconspicuous X-ray signal. 

The deviation of g x from zero by ~ 1 .5<x in the central bin, 
out of the 10 bins we probe, is not unexpected and does thus not 
pose a severe problem for the interpretation of our results with 
respect to cnfw (Sect. l6~4l . 

In a further test, we repeated the analysis centred on Gl, the 
brightest cluster galaxy and found very similar results in terms 
of shapes of g t and g x and fit parameters. 

6.3.2. Likelihood analysis 

While shear profiles serve well to investigate the agreement be- 
tween a cluster shear signal and a mass distribution like NFW, 
there are better methods to infer model parameters, and hence 
the total cluster mass, than fitting techniques. Knowledge of the 
likelihood function 

£=exp (-1 (xtirm, c NF w)-* L (4oo> c^))j (12) 

allows us to quantify the uncertainties in the model parameters 
given the data and - an important advantage over fitting meth- 
ods - also their interdependence. We evaluate the consistency 
between the tangential reduced shear gt,i( r 200>CNFw) predicted 
by an NFW model for the i-th sample galaxy and the tangen- 
tial ellipticity component e t ,, from the data by considering the 
function 



Table 6. Parameters resulting from NFW modelling of 
CL0030+2618 for models relying on different assumptions. 



xl 



z 



|gt,!(/200. CNFw)-£t,i| 

crjU l-|g M (>200,CNFw) 



(13) 



(for a derivation, see ISchneider et ai1 l2000) which we compute 
for a suitable grid of test parameters r2oo and cnfw and determine 
the values r!JS} and c™J^, for which x\ becomes minimal. The 
likelihood approach also allows us to introduce a more accurate 



noise estimator <x; 



SKE 



■ lgt,;l ) f° r the individual shear 
estimate, taking into account the dependence of the noise on the 
shear value itself, as expected by the model PI 

In Fig. |T2j we present the regions corresponding to confi- 
dence intervals of 68.3%, 95.4%, and 99.73% in the r 20 o-c NFW - 
parameter space for three radial ranges in which data galaxies 
are considered. The solid curves denote the fiducial model with 
the complete 0' < 9 < 15' range, giving r™ = 1.52+^gMpc 
and c™^, = 2.0^J j. We will adopt these as the fiducial re- 
sults of our analysis (see Table |6]), yielding a cluster mass of 
Miooir^) = 7.2+2 9 x 10 14 M o (statistical uncertainties) by ap- 
plying Eq. ([8). In the following, we will discuss variations to this 
fiducial case (e.g. the other contours in Fig.fT2l. 

6.4. The Concentration Parameter 

While our resulting r2oo seems reasonable for a galaxy cluster of 
the redshift and X-ray luminosity of CL0030+2618, its concen- 
tration, despite the fact that it is not well constrained by our data 
and cluster weak lensing in general, seems low compared to the 
kno wn properties of galax y clusters: 

iBullock et al] d200ll) established a relation between mass 
and concentration parameter from numerical simulations of dark 
matter haloes, using a functional form from theoretical argu- 
ments: 

cnfw.o / M vir \ 
cnfw = — — | — | (14) 
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5.6+" f° r ^ e v i f i a l mass an d re d- 



1 From Eq. {8}. 

* Mass /J = M2oo/A^20o m umts °f the fiducial cluster mass. 



with cnfw.o = 14.5 + 6.4 and a = -0.15 ± 0.13 for a piv- 
otal mass M, = 1.3 x 10 13 h~ l M Q . This, for z _ 0.50 and 
M vil = 7.2 x 10 14 M© gives c NFW = 3.7. IComerford & Natarajanl 
(120071) . analysing a sample of 62 galaxy clusters for which virial 
masses and concentration parameters have been determined, 
and using the same relation Eq. ( fT4b . find cnfw.o = 9.0 and 
a = -0.13, yielding cnfw 
shiftofCL0030+2618. 

This large interval is consistent within the error bars with 
our fiducial cnfw = 2.0+} 'f with x\l v &of — 13413/13636, but as 
the value itself remains unusually small, we investigate it fur- 
the r. First, we test cnf w = 4.0, close to the value suggested 
bv IBullock eTail (l200lh . while fixing = 1.64Mpc and find 
xtl v Aof = 13423/13637 and the shear profile of the resulting 
model (dash-dotted line in the middle panel of Fig. fTTb to be 
clearly outside the error margin for the innermost bin, demand- 
ing a significantly higher shear in the inner 90" than consistent 
with the measurements. With changes in cnfw mainly affecting 
the modelling of the cluster centre, there is no such tension in 
the other bins. In the next step, we repeat the fit to the pro- 
file, now with cnfw = 4.0 fixed and rioo as the only free pa- 
rameter. The resulting best-fitting model yields = 1.51 Mpc 
(triple-dot dashed in the middle panel of Fig. Qljl, still outside 
but close to the measured lcr-margin of the data. As this fit gives 
AfL,/ v dof - 13418/13637, we conclude that more strongly concen- 
trated models than the fiducial are indeed disfavoured. 

Residual contamination by cluster galaxies reduces the 
measured concentration parameter, as can be seen when "switch- 
ing off" the contamination correction factor (see Table 01. This 
is expected as contamination suppresses the signal most strongly 
in the cluster centre. Removing all galaxies at separations 8 < 
0.5' from the likelihood analysis, we indeed measure a higher 



Use of this error model is denoted writing^ instead of^ 2 . 



c nfw = 2-2+J4, but for the price of larger error bars, as the 
same galaxies close to the cluster centre have the highest con- 
straining power on cnfw- As can be seen from the dashed con- 
tours and the diamond in Fig. Q~2] excising the < 0.5' galax- 
ies just stretches the confidence contours towards higher cnfw, 
leaving r^ 1 = 1 .52+2^ Mpc, and thus the inferred cluster mass 
unchanged (see also Tabled. 

Replacing the contamination correction with a background 
selection down to the faintest magnitudes (jMf a j nt = 29), removing 
a large fraction of the "red sequence-like" galaxies in Fig. [10] 
also yields a higher c?,' = 3.3 ± 1.7 in the shear profile fit, to- 
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Fig. 12. Confidence contours in the NFW pa- 
rameter space spanned by the virial radius r 2 oo 
and concentration cnfw, corresponding to con- 
fidence levels of 99.73%, 95.4%, and 68.3%. 
Also given as symbols are the maximum like- 
lihood cluster parameters for three different 
cases. They are: the fiducial model, using data 
in the complete range 0' < 8 < 15' (solid con- 
tours and dot); a model where the central 30" 



4 g g are excised from the analysis (dashed contours 

and diamond); and a model considering only 
c nfw data inside 0' < 6 < 4' ss r 2 oo (dash-dotted con- 

tours and square). 



• likelihood, fiducial 

♦ likelihood, max( | e | ) = 1 
■ likelihood, max( | e | ) = 1 0"Z 
► likelihood, cen BCG 




Fig. 13. Confidence contours and values of T2oo 
and c NFW maximising the likelihood in depen- 
dence of the maximum shear estimator max(|£|) 
permitted in the catalogue: Given are the fidu- 
cial case (max(|£|) = 0.8, solid contours and dot, 
see Sect. 14. 6t : max(|£|) =1.0 (dashed contours 
and diamond); and max(|fi|) = 10 4 , equivalent to 
no cut (dot-dashed contours and square). Triple 
dot-dashed contours and the triangle denote the 
results for an otherwise fiducial model centred 
on the BCG of CL0030+2618. 



gether with a slightly larger rSL = 1 .67 ±0.18 Mpc and a less sig- 
nificant detection S(0 C ) = 5.10 than the fiducial case. A further 
possible explanation for the low cnfw due to additional lens- 
ing by the z ~ 0.25 foreground structure is rather unlikely (cf. 
Sect. l63~TT l. 

6.5. The Extent of the NFW Profile 

iNavarro et alj (1 19971) designed their profile to represent the mass 
distribution of galaxy clusters in numerical simulations within 
the virial radius. Thus, as theory provides no compelling argu- 
ment to use it out to larger radii, this practice has to be justified 
empirically. 

In the lower panel of Fig. QT| we show results for a toy 
model profile in which the shear signal drops faster than NFW 



outside r2oo- For simplicity, we chose the shear profile of a point 
mass, i.e. 

fext(#)=gt,NFw(#20o)(-^) (15) 

for 9 > 0200, the separation corresponding to r2oo- As in the mid- 
dle panel of Fig. [TT] dashed, dot-dashed, and triple dot dashed 
lines denote the fit to both r2oo and cnfw, setting cnfw = 4.0 
for the same r2oo, and fitting to r2oo for a fixed cnfw = 4.0, 
respectively. The truncation points 6*200 are marked by squares 
in Fig. Q~TJ For the usual two-parameter model with - 
Xl nfw = 0-^0, as for the other two models the truncated, the dif- 
ference in goodness-of-fit between the truncated and pure NFW 
profiles is marginal. 
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Fig. 14. Confidence contours and values of T2oo 
and c NFW maximising the likelihood for mod- 
els of different geometric factor: The fiducial 
case, (Ais/A,) = 0.33 (solid contours and dot) 
is compared to the low ({D^/D,,} = 0.30, 



dashed contours and square) and high cases 
((Ais/A) = 0.36, dot-dashed contours and di- 
amond) derived from the CFHTLS Deep fields 
in Sect. lB32l 



Secondly, we repeat the likelihood analysis for galaxies 
0' < < 4' =s 0200 only. The dash-dotted contours and the square 
in Fig. [12] for the resulting optimal parameters show the corre- 
sponding values. Here, = 1.58+'^ Mpc and c™ n w = 1.6+^ 
are more degenerate than in the fiducial case (cf. Table [6]). We 
conclude that there is no evidence in the CL0030+2618 data for 
a deviation of the shear profile from NFW at r > 7-200. Applying 
Occam's razor, we use this profile for the whole radial range, but 
stress cautiously that we cannot preclude an underestimation of 
the errors and, to a lesser extent, a bias in the virial mass here. 

6.6. Shear Calibration 

As already pointed out in Sect. 14.61 the maximum shear esti- 
mator max(|e|) considered in the catalogue strongly affects aver- 
aged shear observables. In Fig. [13] we quantify this dependence 
by comparing the confidence contours and best values for r2oo 
and cnfw from the fiducial max(|e|) = 0.8 catalogue (solid con- 
tours and dot) to cases with max(|e|) = 1.0 (dashed contours and 
diamond) and max(|e|) = 10 4 (dot-dashed contours and square). 
The latter includes even the most extreme shear estimates^. The 
max(|e|) cut, via the amplitude of the shear signal, mainly influ- 
ences r^j ( ", reducing it by 6 % (13 %) for the frequently used 
max(|e|) = 1.0 and the extreme max(|e|) = 10 4 , respectively. In 
turn, the mass estimate would be reduced by 17 % (35 %), as 
can be seen from Table [6] 

The influence on the mass estimate by the choice of 
max(|e|) is compensated by the shear calibration fo + 1 and one of 
the effects which we account by considering different fo. Given 
the uncertainty <x fo = 0.05 (Sect. 14.61 ), we repeat the likelihood 
analysis with fo = 1.13. For the negative sign, the signal dilution 
by foreground galaxies has to be taken into account. Combining 
in quadrature the 1 8 % foreground dilution estimated from the 
CFHTLS Dl field (Sect. lB~XTT i with cr fo , we arrive at f = 0.88 



12 Note that, although unphysical, shear estimates £> 1 in KSB are to 
some extent justified when averaging over large ensembles. 

13 The sign here is likely due to a statistical fluke; theory expects rSJJ 
to increase with a less strict max(|e|). 



as the lower bound of the error margin. The +5 % (-20 %) varia- 
tion in fo translates into +1.3 % (-7.2 %) in r^, yielding again 
» +4 % (» -20 %) variation in M 20 o (see Table©. 

6. 7. Combined Mass Error Budget 

Replacing the weak lensing centre in our fiducial model by the 
cluster's BCG as the centre of the NFW profile, we find the re- 
sulting differences in r^j ( " and c^p W returned by the likelihood 
method, and hence in M200, to be small (cf. triple dot-dashed 
contours and triangle in FigfO] Table©. We conclude the error 
on the chosen centre to be subdominant. 

Variations in the geometric factor {D ds /D s ) induce a similar 
scaling in rSjj and c™ w as shear calibration does. Using the er- 
ror margin from the determination of the distance ratios from the 
CFHTLS Deep fields (Sect. IB.3.2t . we produce likelihood con- 
tours for (Dds/-D s ) = 0.30 (dashed lines and square in Fig.PBli and 
(Z)ds/£> s ) = 0.36 (dot-dashed contours and diamond). Comparing 
to the fiducial model (solid contours and dot), we find an in- 
crease in rl™ by 4.6 % and by 14 % in M 20 o for (D d JD s ) = 0.30 
(a more massive lens is needed for the same shear if the source 
galaxies are closer on average) and a decrease by 4.6 % in 
and 13 % in M 200 for (D d JD s ) = 0.36 (cf. Table©. 

An additional source of uncertainty in the mass estimate 
not discussed so far are triaxiality of galaxy cluster dark mat- 
ter haloes and projection of t he lar ge-sc ale structure (LSS) onto 
the image. King & Corlessl (120071) and [Corless & Kingl (l2007h 
showed with simulated clusters that masses of prolate haloes 
tend to get their masses overestimated in weak lensing while 
masses of oblate haloes are underestimated. 

Again owing to cosmological simulations. iKasun & Evrardl 
(2005) devised a fitting formula for the largest-to-smallest axis 
ratio 77 of a triaxial haloes as a function of redshift and mass 

.(M2oo,z)^o(l + zr(l-^ln(^^)) (16) 

with e = 0.086, £ = 0.023, and 770 = 0. 633. I nserting the values for 
CL0030+2618, we find 77 = 0.61 and, like iDietrich et all d2009h 
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Table 7. Parameters and fit values of the |Vikhlinin et al. (2006) 
ICM model for CL0030+2618. 



Parameter 


Quantity 


Fit Value 


"o 


pivot density 


3.784 x 10- J cnr J 


r c 


core radius 


139 kpc 


r s 


scale radius 


421 kpc 


a 


exponent 


0.5867 


P 


exponent 


0.4653 


E 


exponent 


1.2293 



whose lin es we are following, der ive the following maximal bi- 
ases from ICorless & King] (|2007): for a complete alignment of 
the major cluster axis with the line of sight mass is overesti- 
mated by 16 %, while complete alignment with the minor axis 
results in a 10 % underestimation. 

The projection of physically unrelated large scale structure 
can lead to a signifi cant underestimation o f the statistical errors 
in M?on and cn fw dHoekstral I2003L 120071) . The simulations of 
iHoekstral d2003l) yield an additional error of + 1 .2 /T I xlO I4 M = 
± 1 .67x 1O 14 M for a cluster in the mass range of CL0030+26 18, 
and little redshift dependence for z > 0.2. Thus, we adopt this 
value as the systematic uncertainty due to large scale structure. 

We define the systematice mass uncertainty cr sys as the 
quadratic sum of the errors o- ca \i from shear calibration, cr geom 
from the geometric factor, cr pl0 j from projection, and <x L ss from 
large-scale structure^ The total error, used in Fig. Q3] is then 
defined as the quadratic sum also including <x stat : 



sys 



2 , 2 



LSS" 



pro] 



cali 



(17) 



We note that the statistical errors are already quite large and 
the dominating factor in Eq. {T7J. As its main result, this study 
arrives at a mass estimate of M 2 oo = l-l-lt-l'l x 10 14 M© for 
CL0030+2618, quoting separately the statistical and systematic 
error as the first and second uncertainty. 

6.8. Comparison to X-ray masses 

We will now compare this weak lensing mass to a mass profile 
drawn from the Chandra analysis of CL0030+2618. Under the 
assumption that the ICM is in hydrostatic equilibrium, the total 
mass M(< r) of a ga laxy cluster within a radius r can be derived 
as fcf. lSarazirIT9 88): 



M(;r)= -W)r ( dhp g | dln7 x 



pm ? G 



d In r d In r 



(18) 



where m p is the proton mass and /i the mean molecular mass. In a 
first step, we treat the ICM temperature to be independent of the 
radius and fix it to the lVikhhnin et al.l d2009al) value of k B (Tx) = 
5.63 + 1.13 keV\ For the gas density p g , we use a (simplified) 
IVikh linin et al. (2006) particle density profile 



V«e«p 



(r/r c T 



(l +r 2/ r f- a/2 (l +r r/rj) 



(19) 



with the parameters given in Table [7] and a fixed y = 3 and arrive 
at a mass of M hy d = (6.26 ± 1 .26) x 10 14 Mq at the virial radius 
of f*2oo = 1-52 Mpc obtained in the lensing analysis. We show the 



14 We remark, however, that strictly speaking o" LS s qualifies as a sta- 
tistical error. 



corresponding mass profile as the thick and its error margin as 
the grey lines in the upper panel of Fig. Q3] 

This value is in very good agreement with the weak lensing 
mass estimate (dot with thick error bars for statistical and thin 
error bars for systematic plus statistical uncertainties in Fig. [T3Tl. 
The consistency between the X-ray mass profile derived from Tx 
and the (baryonic) ICM using Eq. [19] and the NFW profile de- 
scribing the combined dark and luminous matter densities holds 
at all relevant radii > 50 kpc in a wide range from the cluster core 
till beyond the virial radius. 

Assuming an isothermal cluster profile, one likely overes- 
timates the total hydrostatic mass, as the ICM temperature is 
lower at the large radii dominating the mass estimation around 
r v ; r . The competing effect of the temperature gradient term in 
the hydrostatic equation is subdominant compared to this effect 
of the temperature value. 

Therefore, to estimate the systematic uncertainty arising 
from isothermality, we consider a toy model temperature pro- 
file consisting of the flat core at (Tx), an power-law decrease 
at larger radii, and a minimal temperature k B To = 0.5 keV in 
the cluster outskirts to qualitatively represent the features of an 
ensemble-averaged temperature profile: 



k B T x (r) = 



k B (T x ) r< n 
p r 9 rj<r<r t 
k B T r>r x 



(20) 



where we choose a core radius = r2oo/8 (as used in lPratt et al.l 
2007), a power-law slope q = -0.4 taken as a typical value 
found by Eckmiller et al. (in prep.), and fixing the truncation 
radius r t and amplitude p demanding continuity of Tx(f). The 
mass profile resulting from this temperature distribution is plot- 
ted in Fig. [15] (upper panel) as the dash-dotted line, giving an 
estimate of the systematic uncertainty in the X-ray profile. Its 
value coincides with the lower end of the lcr mass range for M w i 
at r2oo, taking into account its systematic errors. Another sys- 
tematic factor in X-ray analysis is non-thermal pressure support, 
leading to an unde restimation of the X-ray mass by ~ 10 % (e.g. 
IZhang et al.l2008l) . Taking into account all these effects, we con- 
clude a very good agreement of X-ray and weak lensing mass es- 
timates of CL0030+2618, despite the potential perturbation by 
the line-of-sight structure. 

In the lower panel of Fig. [15] we show the ratio Mh y d/M w i 
of hydrostatic X-ray and weak lensing mass as a function of ra- 
dius. Although this quantity has a large error, our values are in 
good agreement w ith the X-ray-to-lensing mass ratios found by 
IZhang et al.l(l2010h for a sample of relaxed clusters for three radii 
corresponding to overdensities A = Pnfw(< r)/p c = 2500, 1000, 
and 500 (black line). We note that we recover well the relation 
^f(A) found bv lZhang et all (l201dh by fitting their cluster sam- 
ple data (grey line). 



7. Summary and Conclusion 

With this study, we report the first results for the largest weak 
lensing survey of X-ray selected , high-redshift c l usters , the 400d 
cosmological sample defined bv lVikhlinin et al.ld2009al) and de- 
termine a weak lensing mass for an interesting cluster of galax- 
ies, CL0030+2618, which had not been studied with deep opti- 
cal observations before. We observed CL0030+2618, along with 
other clusters of our sample, using the Megacam ~ 24' x 24' im- 
ager at the MMT, ob taining deep g'r'i' exposures. Employing 
an adaptation of the lErben et al] d2005h pipeline, THELI, and 
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Fig. 15. Comparison of mass profiles of 
CL0030+2618. 

Upper panel: The hydrostatic mass M hyd (< r) 
derived from the Chandra analysis (thick solid 
line). A constant ICM temperature is assumed 
and the grey lines delineate the error margin de- 
rived from its error. The dash-dotted line gives 
the Chandra profile for a more realistic tem- 
perature profile. The dot with error bars and 
the dashed line denote the mass estimate and 
profile Af w i(< r) from our weak lensing analy- 
sis, assuming an NFW profile. The thick error 
bars show the statistical errors while thin bars 
include all components dicussed in Sect. 16.71 
Lower panel: Ratio of X-ray to lensing mass 
as a function of radius (black line). The sym- 
bo ls and grey line sh ow the Mh y i/M„\ found 
by Zha ng et al] d2010h at three overdensity radii 
and their fitted relation. 



the "TS" KS B shape measuremen t pipeline presented by T. 
Schrabback in Hev mans et al.l d2006l) . we, for the first time, mea- 
sure weak gravitational shear with Megacam, showing its PSF 
properties to be well suited for such venture. 

The lensing catalogue of background galaxies is selected by 
a photometric method, using g'r'i' colour information. Despite 
similar number count statistics, we find different photometric 
properties in our Megacam field than in the CFHTLS Deep fields 
used to estimate the redshift distribution of the lensed galaxies. 
The photometr ic measurements establish the galaxy we name 
Gl, for which lBovle et al.l ( 1 19971) determined a redshift z- 0.5 16 
as the BCG of CL0030+2618, ruling out a slightly brighter 
source found inconsistent in its colours with the cluster redshift 
z = 0.50. We find additional evidence for the presence of a fore- 
ground structure at z ~ 0.25 from photometry but find it does 
neither significantly affect the lensing nor the X-ray mass esti- 
mate of CL0030+26 18. 

Having applied several consistency checks to the lensing 
catalogue and optimising the aperture mass map of the cluster, 
we detect CL0030+2618 at 5.84cr significance. The weak lens- 
ing centre obtained by bootstrapping this map is in good agree- 
ment with the BCG position and the X-ray detections by Rosat, 
Chandra, and XMM-Newton. Two tentative strong lensing arcs 
are detected in CL0030+2618. 

Tangential alignment of galactic ellipticities is found to 
extend out to 10' separation and well modelled by an NFW 
profile out to > 2r2oo- The low concentration parameter found 
by least-squares fitting to the shear profile is confirmed by the 
likelihood method with which we determine CL0030+2618 to 
be parametrised by r™ = 1 .52+°;g Mpc and c™ w = 2.0^. 
Modifying the likelihood analysis for the fiducial case, we es- 
timate systematic errors due to shear calibration, the redshift 
distribution of the background galaxies, and the likely non- 
sphericity of the cluster. Further, we confirm the best model to 
change little if the BCG is chosen as cluster centre rather than 
the weak lensing centre. We arrive at a virial weak lensing mass 
for CL0030+2618 with statistical and systematic uncertainties 
of M 2 oo = 7.2+2 9-2 5 x 10' 4m O> in excellent agreement with the 



virial mass obtained using Chandra and the hydrostatic equa- 
tion, M hyd = 6.26 + 1 .26 x 1O 14 M . 

Noting that statistical errors on the lensing mass are still 
high, we conclude that high-quality data and well-calibrated 
analysis techniques are essential to exploit the full available cos- 
mological information from the mass function of galaxy clusters 
with weak lensing. Nevertheless, once lensing masses for all the 
36 clusters in the sample are available, these statistical errors are 
going to be averaged out and reduced by a factor of 6 by combin- 
ing all clusters when testing for cosmological parameters. Thus, 
understanding and controlling systematic errors remain impor- 
tant issues. We are going to proceed our survey with the analy- 
sis of further high-redshift clusters from the 400d cosmological 
sample observed with Megacam. 
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Appendix A: Details of Data Reduction 

In this Appendix, we provide a detailed view of the subtleties of 
the data reduction outlined in Sect. [3] 

A.1. Chips and Amplifiers 

The MMT Megacam control software offers a number of options 
for the CCD readout. As already mentioned in Sect. l2.2l there are 
36 physical CCD chips each of them is equipped with two output 
amplifiers, giving a readout of 1024 x 4608 (unbinned) pixels 
per amplifier. For our programme, we have chosen to use all 72 
amplifiers each reading out half a chip, thus reducing readout 
time by a factor of two. As a result, Megacam raw images are 
multi-extension fits files with 72 extensions. 

Owing to this, all run processing tasks are performed on 
the 72 subframes individually. Files from the two chips of an 
amplifier are joined at the end of the run processing prior to the 
astrometric calibration in order to increase the usable surface for 
the astrometric procedures. 

A.2. The "Run Processing" Stage 

- De-biasing: By stacking all bias frames taken within a suit- 
able time interval around the date of science observations, 
a master bias image is constructed and subtracted from all 
other frames. 

- Flatfielding: THELI applies a two-step process. First, sci- 
ence frames are divided by a master sky flatfield frame. In 
the second step, the median of all science frames is calcu- 
lated, discarding the pos itions at which object s have been 
detected by SExtmctor dBertin & Arnoutsl[T996l) . Due to the 
dithering, for every pixel in the field-of-view, these "super- 
flats" contain signal from the sky background from slightly 
different positions on the sky. Thus, the superflat provides a 
measure to compare the response in different pixels. 
Selecting the frames to contribute into the superflat to 
achieve the optimal flatness of the background is the most 
time-consuming and work-intensive step in run processing, 
as inhomogeneities in individual frames will propagate into 
the superflat. Imperfect photometric conditions and vari- 
able instrumental gains are two common reasons for science 
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frames to be removed from the calculation of the superflat. 
Very bright stars near target clusters, e.g. HIP 9272 (BD- 
00301) with V = 8.28 at 2' distance from CL0159+0030, 
exacerbate the situation. Involving many iterations of the - 
manual - frame selection process, our superflatfielding is ef- 
fective in reducing the relative background variation over 
the field in the superflatfielded exposures to < 1.5 %, and 
to < 1.0 % for most exposures. 

In the superflatfielding stage, the different sensitivities of the 
amplifiers are determined and equalised, taking into account 
all exposures within the THELI run. This we can do, be- 
cause the relative sensitivities of most of the amplifiers are 
constant most of the time. Gain equalisation is achieved by 
scaling each ampl ifier with an appropriate factor detailed in 
lErben et aU d2005. Sect. 4.7). Some amplifiers, however, ex- 
perience gain fluctuation on short timescales of the order of 
days. In such situation the same superflatfield frame can no 
longer provide the same quality of flattening to all exposures 
such that we had to split the g'-band data taken on Nov. 8, 
2005 from the rest of the exposures taken on Oct. 30, 2005 
and Nov. 1, 2005, at the cost of a lower number of exposures 
contributing to the superflatfield in each of the two sub-runs. 

- De-fringing: In the bands where this is necessary, the fring- 
ing pattern can be isolated from the high spatial frequen- 
cies of the superflat and subtracted from the science frames. 
On the other hand, we divide by the superflat containing the 
lower spatial frequencies which carry information about the 
(multiplicative) "flatfield" effects. 

- Satellite tracks: We identify satellite tracks by visual in- 
spection when assessing frames for superflat construction 
and mask pixels which are affected in the given exposure. 
Masked pixels (stored as a DS9 region files) are set to zero in 
the construction of the weight images. 

- Weight images: Taking into account bad pixel information 
from the biao, flatfield, and superflatfield frames we con- 
struct a weight image, i.e. noise map, for each individual 
amplifier and exposure in the run. Our algorithm is not only 
sensitive to cold and hot pixels but also to charge "bleeding" 
in the vicinity of grossly overexposed stars. 

A3. The "Set Processing" Stage 

- Astrometry: We perform the astrometric calibration of 
our data using the best available catalogue as a refer- 
ence. In case of overlap with SDSS Data Release Six 
dAdelman-McCarfhv et"al]|2008l) this is the SDSS catalogue; 
otherwise we emp loy the shallower USN0 Bl catalogue 
dMonet et al.l 12003). as it is the densest all sky astromet- 
ric catalogue. The astrometric calibration is carried out by 
the TERAPIX software Scamp dBertinll2006l) . replacing the 
Astrometrix programme earlier used within THELI. We find 
Scamp to be more robust than Astrometrix when working 
on chips with a small field-of-view on the sky, as given for 
Megacam (see Sect. [3]). Compared to the otherwise simi- 
lar design of MegaPrime/MegaCam at the Canada-France- 
Hawaii Telescope the MMT Megacam chips cover ~ 1/6 of 
the solid angle on the sky, reducing the number of usable 
sources for astrometry by a similar factor leading to less ac- 
curate and robust astrometric solutions in case these are cal- 
culated on a chip-to-chip basis. 



IS Here, we also use dark frames, although they are not necessary for 
running THELI. 



The most important innovation is, that while Astrometrix de- 
termines an astrometric solution for each chip (Megacam am- 
plifier in our case) individually, Scamp recognises the am- 
plifiers of one exposure belonging together and can take into 
account information on the array configuration, drastically 
reducing the effort to be invested into this task. 
We provide these additional constraints by defining a tem- 
plate for the same instrument configuration and filter. This 
template is drawn from the observation of a dense field, i.e. 
a star cluster. This template guarantees a sensible solution 
even with few (< 20) astrometric standard stars per chip, a 
condition frequently met with Megacam in poor fields. 
Furthermore, by running Scamp on all frames in all filters 
for a given target cluster with only one call to the software, 
we ensure consistency of the astrometric solutions among the 
THELI sets corresponding to the resulting stacks in different 
passbands. 

For the combined data set of CL0030+2618, we achieve an 
accurate calibration with a lcr intrinsic accuracy of 0704 of 
the sources detected with Megacam and 0727 with respect to 
the astrometric standard catalogue USN0 Bl. 

- Relative photometry: Simultaneously with the astrometric 
calibration, the relative photometric zeropoints of the frames 
are established by Scamp. In the first part of this two-step 
process, relative zeropoints are determined only from the dif- 
ferences in flux found for the astrometric reference stars in 
different exposures. These are independent from the absolute 
photometric calibration detailed in Sect. 13.21 

In this first step, fluxes of the same object in different expo- 
sures are compared. For the coadded image resulting from 
stacking to be well-calibrated, the variation in relative zero- 
points among the contributing frames needs to be small. We 
decide to include only images which show a zeropoint less 
than 0. 1 magnitudes from the median zeropoint: 

Z re i-median(Z re i)<-0.1 (A.l) 

In the second step, if the absolute photometric calibration 
(Sect. 13 -2b has been applied already, we comp u te the cor- 
rected zeropoints defined in Hil debrandt et al.l (120061 Eq. 
(2)) of those individual frames we consid er to be taken un- 
der ph otometric conditions. As detailed in Hilde brandt et al.l 
(2006), corrected zeropoints are a useful consistency check, 
as they are the same for exposures obtained in photometric 
conditions. 

- Coaddition: Conforming with THELI standard, SWarp is 
used to stack ("coadd") images. This, together with the 
Scamp astrometry, also removes optical distortions, yield- 
ing a constant pixel scale in the coadded image. The final 
products of the set stage are the coadded image and the cor- 
responding weight image (Fig. [2]). 

A.4. Image Selection 

The success of a lensing analysis crucially depends on the data 
quality. Because of the necessity to establish a common image 
coordinate system and to rebin all data onto the new grid, the 
stackingprocess is a potential source of biases to the shape infor- 
mation[3 It is evident that the decision which frames should con- 
tribute to the shape measurement is of great importance. Apart 
from seeing and photometric quality which can be easily as- 



16 There is an ongoing debate whether shapes should be measured on 
individual frames, instead. 
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Fig. A.2. Spatial distribution of stellar anisotropics for an example exposure of high overall PSF anisotropy. Shown are the sizes 
and orientations of the raw ellipticity e for stars identified in the MMT Megacam exposures of CL0030+2618 labelled Q936 {upper 
panel) and 0952 (lower panel) in Fig. I A. 1 1 While within each chip the x and y pixel axes are to scale the array layout is only 
schematic. 
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Table A.l. Coefficients of photometric calibration defined by Eq. © for the photometric nights used to calibrate the observations 
ofCL0030+2618,CL0159+0030,andCL0809+2811. 



Filter 


Obs. Date 


: T 




Pi 


colour index 


Yt 


: T 

"par 


g' 


2005-10-30 


27.27 ± 


0.02 


0.106 ± 


0.007 


g '-r> 


-0.14 + 0.02 


3 






97 1 S + 


n n? 


0.124 ± 


0.008 


g ~ r 


— n ok + n ni 


T. 

J 




2005-11-01 


27.35 ± 


0.02 


0.115 ± 


0.005 


g'-r' 


-0.21 ± 0.02 


3 


i' 


2005-10-30 


26.49 ± 


0.02 


0.127 + 


0.004 


r' — i' 


-0.11 ± 0.02 


3 




2005-10-31 


26.47 ± 


0.01 


0.122 + 


0.002 


f - i' 


-0.09 ± 0.01 


3 




2005-11-01 


27.41 ± 


0.01 


0.119 + 


0.002 


r' - V 


-0.03 ± 0.01 


3 


r 1 


2005-10-30 


26.95 ± 


0.02 


0.046 ± 


0.002 


g'-i' 


-0.10+0.02 


3 




2005-10-31 


26.90 ± 


0.01 


0.043 ± 


0.003 


g'-i' 


-0.05 ± 0.01 


3 




2005-11-01 


26.96 ± 


0.01 


0.048 ± 


0.004 


g'-i' 


(-0.10) § 


2 




2005-11-08 


26.81 ± 


0.01 


0.046 ± 


0.003 


g'->' 


(-0.10) § 


2 



' Normalised to an exposure time of Is and an airmass a = 0. 
* Number of parameters used in the fit. 
§ Fixed to default value. 



sessed while the observation takes place, PSF anisotropy is a 
key factor as it can only be corrected up to a certain degree. 

We found that irrespective of the seeing, a significant 
fraction of the Megacam exposures from the observing runs 
in October 2004, October 2005, and January 2008, covering 
the four clusters CL0030+2618, CL0159+0030, CL0230+1836, 
and CL0809+281 1 suffer from highly elliptic PSF which shows 
little variation over the field-of-view. 



Inspecting the anisotropy stick plots like Fig. lA.2l for a fair 
fraction of all frames taken in these three above-mentioned runs, 
we come up with the following criterion: If the mean ellipticity 
in the stars is (|e|)<0.05 (the smaller circle in Fig. IA.lb the vari- 
ations over the field due to the properties of the optical system 
are clearly discernible. These frames do not suffer from a large 
tracking error and thus we include them in the analysis in any 
case. On the other hand, exposures whose stars are more ellip- 
tical than (|e|) > 0.06 (the larger circle in Fig. I A. 1 b show a stel- 
lar anisotropy which is mostly constant over the field which we 
can attribute to MMT tracking errors of varying strength. These 
frames are excluded from the lensing analysis. In the intermedi- 
ary case of 0.05 < (|e|) < 0.06 we decide on a case-to-case basis 
by inspecting the respective anisotropy stick plot where frames 
in which the "tracking error-like" contributions seem prevalent 
are excluded. 



A.5. The Megacam filter system 

To establish the transformation between MMT and SDSS mea- 
surements, we need to know the transmissivities of both instru- 
ments to a great detail. For Megacam, the instrument websiteTI 
offers detailed laboratory transmission curves of the actual fil- 
ters and a few data points on the CCD's quantum efficiency. 
We average the tabulated quantum efficiency values over the 36 
Megacam chips. Correspondingly, the SDSS collaboration pro- 
vides data on the combined sensitivity of its camera/filter sys- 
terrQ Assuming the atmospheric absorption to behave the same 
on both sites, we can directly compare the responses of the two 
instruments, as visualised in Fig. [A3] 



17 Overview: 

http : //www . cf a . harvard . edu/mmti/megacam . html; 
filter data: 

http : / /www . cf a . harvard . edu/"bmcleod/Megacam/Filter s/ 

18 http: //www. sdss.org/dr7/instruments/imager/ 



A.6. Results of Photometric Calibration 

Photometric calibration is achieved by fitting Eq. © to 
the instrumental magnitudes of the photometric standards 
(Sect. [3727Tb . For each filter, we chose a colour index in Eq. <j4j 
which has been proven a reliable transformation in calibra- 
tion of the Canada-France-Hawaii Telescope Legacy Survey 
(CFHTLS) data which also uses a similar filter systerro These 
colour indices are given in Table IA.3I which shows the results 
for the fit parameters Zf, /3f, and jf for the photometric nights 
used to calibrate the CL0030+2618 data (i.e. the datasets for 
CL0030+2618, CL0159+0030, and CL0809+2811; Sect.[372"2}. 

Comparing the colour terms pf for the different nights, we 
find considerable agreement of the values for all three bands, al- 
though the error bars from fitting Eq. (HJl might underestimate 
the true errors. While our determinations of fi r are all consistent 
with each other, there is som e tension between th e f3 g values. In 
previous Megacam studies, (lHartman et al.l l2008. Table 5) quote 
(3 g = 0.122 ± 0.002 and # = 0.137 + 0.002, the first in agree- 
ment to our r esults, the lat ter significantly higher than our value. 



Furthermore, Wals h et alJ d2008) find yS g =0.091 + 0.068, consis 



tent with our values due to the large error bar. We suggest that 
the large span in values for /3 g might be due to the dependence of 
the filter throughput as a function of distance to the optical axis, 
which is strongest an this band. Further investigation is needed 
to conclude about this issue. 

Appendix B: The Background Catalogue 

BA. Photometric Analysis: The Red Sequence 

In clusters of galaxies at low and moderate redshifts, early-type 
galaxies, i.e. elliptical and spheroidal systems tend to be more 
common than disk galaxies. Cluster galaxies are observed to 
be deficient in gas and thus show little ongoing star forma- 
tion but a re on average dom inated by old, red stellar popula- 
tions (e.g jBower et al.lll992h . In fact, cluster galaxies represent 
the reddest galaxies observed at a given redshift with what are 
considered the most gas-depleted system s showing very simila r 
colours over a large range in magnitude (Gladd ers et alJ [l998). 
Observationally, this cluster red sequence is one of the currently 
most prolific m ethods in detecting clusters o f galaxies in the op- 
tical band (e.g. JGladders & Yeell2000ll2005h . 

We consider the (g'-i') vs. i' colour-magnitude diagramme 
of the galaxies in the galaxy shape catalogue (i.e. before cuts to 



http://www3.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/megapipe/docs/filters.html ) 
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Fig. B.2. Colour-colour selection of the lensing catalogue: plotted are the g'—f vs. r'-i' colours of the objects in the galaxy shape 
catalogue (with cuts on min(trP g ), min(v ga i), and max(|(e|) already applied). The galaxy sample is divided into three magnitude 
bins by the m^ght an d mfaint parameters. All sources brighter than ;«bdght (largest dots) are removed in producing the final lensing 
catalogue while all sources fainter than Wf a ; nt smallest dots) are kept. Intermediately bright galaxies with m^ght <r' < mfaint (medium- 
sized dots) mark the transition between these two regimes. Only in this magnitude interval, the selection into the final galaxy 
catalogue by colour indices applies: sources outside the thick polygon bounding the region in which we find the brighter galaxies 
and likely cluster members are included in the final catalogue. See Table IB. fl for the definition of the polygon tracing the locus of 
bright galaxies. 



Table B.l. Tested regions in colour-colour space inside which galaxies with m^ght < r ' < '«faint are excised from the lensing 
catalogue. 



Polygon 


min (r'-i') 


max (r'-i') 


min (g'-r') 


max (g'-r') 


min (s) 


max (s) 


large 





1.25 


0.25 


2.25 


-1.0 


0.583 


small 


0.1 


1.2 


0.4 


2.0 


-0.733 


0.1 



select sources for their lensing signal apply) close to the coor- 
dinates of CL0030+2618 in order to identify the red sequence 
of this z = 0.5 cluster as the observed g' and i' passbands are on 
different sides of the Balmer break at the cluster redshift. 

Having removed most extended galaxies early-on in the 
KSB pipeline, we do not expect to find the most prominent clus- 
ter members in the catalogue for which shear estimates are deter- 
mined. Indeed, the upper panel of Fig. lB.ll shows a rather broad 
distribution in g' - i' colour of the galaxies at 9 < 4' from the 
Rosat cluster centre. Nevertheless, we find an enhancement in 
the number of galaxies extending from around (g'-i') ~ 2.8 for 
the brighter (/' w 21) to (g'-i') ~ 2.3 for the fainter (z" * 27) 
sources in our catalogue, especially coming from a high number 
of galaxies very close (6<2') to the cluster centre. 



The CWW80 template for an elliptical z = 0.50 galaxy pre- 
dicts (g'-i') ~2.8. This (solid line and large dot at i' = 20 in the 
upper panel of Fig. lBTTT i is in good agreement with the bright end 
of our observed tentative cluster red sequence, indicating we in- 
deed detect the red sequence of CL0030+2618E3 At z = 0.25, the 
tentative redshift of the foreground structure, the same template 
yields (g'-i') ~ 2.1 (dashed line and large dot at i' = 20 in the 
upper panel of Fig. IB. II ). The broad distributions in g'-i' colours 
and the not very distinctive red sequence of CL0030+2618 are 
in agreement with the presence of a foreground group. 



20 For this argument, we can neglect the known slope of the red se- 
quence due to the lower me tallicity of the many d warf galaxies among 
the fainter cluster members (Glad ders et al.| [l998). 
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Fig. B.3. Upper panel: The maximum M ap signal-to-noise ratio found in the vicinity of CL0030+2618 as a function of the back- 
ground selection introduced by mf a i nt and mtnght- The solid line corresponds to a magnitude cut m^ght = m faint while the dotted, 
dashed, dot-dashed, triple dot-dashed, and long-dashed lines denote background selections by galaxy colours in increasingly wide 
intervals of mf a i nt -mbri g ht = {0.5, 1.0, 1.5, 2.0, 2.5} respectively. Here, the smaller polygon in Fig. IB.2l is used, assuming a well-defined 
locus of cluster galaxies in colour-colour-space and, in turn, for a rather inclusive selection of galaxies. Lower panel: The number 
N of galaxies in the shear catalogue as a function of mbright and mf a ; nt . The horizontal line gives N ~ 16000 before applying any 
background selection for comparison. The colours and line-styles denote the same cases as in the upper panel. 



In the lower panel of Fig. IB. II we show the g' - r' colours 
of the same central region galaxies as a function of their r' - i' 
colours. (Compare also Fig. IB. 21 ) In addition to the main clump, 
there is a population of galaxies with both red g' - r' and r' - i' 
colours following the locus of the bright galaxies in Fig. IB. 21 
As expected, the CWW80 templates for z ~ 0.25 and z ~ 0.50 
belong to the redder population, and would for m < mf a j nt be 
excluded from the lensing catalogue by the background selection 
(Sect.lB~2l. 

B.2. Background Selection by Galaxy Colours 

The selection of background galaxies is based on r'—i' vs. g'—r' 
colour-colour-diagrammes for galaxies of intermediate magni- 
tude works as follows: We identify the region in the colour- 
colour-diagramme populated by the brightest galaxies, a sample 
we assume to be dominated by cluster ellipticals. As the clus- 
ter red sequence shows, the colour of early-type systems in a 
cluster of galaxies varies only little with magnitude. Indeed, as 
can be seen from Fig. IB. 21 the bright galaxies observed in the 
CL0030+2618 field show a well-defined relation between their 
r'—i' and g'—r' colours. Inferring that the fainter cluster members 
which fall into the mbright < r' < mf a ; nt interval on average show 



the same colours as their brighter companions, we exclude those 
intermediately bright galaxies with colours similar to those of 
the brighter objects while keeping those inconsistent with the 
colours of the bright sample. 

Following a me thod introduced bv lBradac et al. (2005]) and 
iKausch et alj {2007), we empirically define two polygonal re- 
gions in the colour-colour diagramme, a "small", rather inclusive 
polygon and a "large" polygon for a more conservative selection 
(thick and thin lines in Fig. IB. 21 respectively). We test the in- 
fluence of the colour-colour selection on the lensing signal for 
those two cases. Table |B~T1 gives the respective limits in g' — r 1 , 
r'—i' , and the second-order colour index s :- %r'—%i'—g' chosen 
to be parallel to the locus of the bright galaxies in Fig. IB. 21 

Figure IB. 31 (upper panel) shows the effect of the back- 
ground galaxy selection on the S -statistics if the "small" poly- 
gon defined in Table IB. II is used for the intermediary bright 
galaxies. Here, the solid line denotes a pure magnitude cut at 
'"bright = '«faint while the different line styles show cases in 
which the colour-colour criterion acts in different intervals of 
'Wfaint - mbright- First thing to note is that the S -statistics depends 
more sensitively on mf a ; nt than on 'weight, w i m i ts maximum in 
the range 22.0 < mf a i nt < 22.5, irrespective of mbright- The greater 
relative importance of m^ mt does not come as a surprise as, in 
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the r' < 25 magnitude range we study here, source counts are 
rising steeply towards fainter magnitudes (Fig. IB.4t . 

Secondly, we notice that the improvement in the S -statistics 
upon using the best value of mf a ; nt = 22.5, which we now adopt, 
over the case of not applying photometric criteria to our cata- 
logue (corresponding to mf a i nt = 17.6) is small: S = 5.73 for 
'"bright = '"faint as compared to S = 5.46. This may partly be ex- 
plained by the small number of catalogue objects affected by 
background selection. As can be seen by comparing the num- 
ber of objects in the lensing catalogue as a function of mf a i nt 
and mbright in the lower panel of Fig. IB. 3 1 with the S -statistics, 
as selection starts removing (signal-diluting foreground) galax- 
ies from the catalogue at mf ault > 21.5, the S -statistics begins to 
increase around the same point. For instance, with a magnitude 
cut at ;«f a i n t = 22.5, the remaining 92.5 % of the sources yield a 
S = 5.73, while for a mf a j nt = 21 .5 magnitude cut, the remaining 
97.3 % of the catalogue give S =5.53. 

Similarly, the strong decrease in detection significance for 
»if a int > 22.7 - most pronounced for the '«bri g ht = "ifaint case - can 
be explained by a cut at faint magnitudes rejecting an increas- 
ingly large number of signal-carrying background galaxies. For 
the various 'Wbright < '"faint cases, the higher signals for a given 
'Mfaint demonstrate that galaxies of intermediary magnitude with 
colours inconsistent with cluster ellipticals are kept in the cata- 
logue and contribute to the signal. 

Repeating this analysis with the "large" polygon defined in 
Table IB. II we find the dependence of S on '^bright for a given 
Mfaint to be largely reduced. This can be explained by the re- 
strictive choice of the "large" compared to the "small" polygon, 
leaving only few galaxies of intermediary magnitude in the cat- 
alogue. 

For the following analyses, we choose the "small" polygon 
and the parameter combination mf a i nt = 22.5, ntbright — 20.0, yield- 
ing the near-optima^]] overall detection of the cluster: S = 5.84. 
We also tested catalogues with 'Mfaint - '"bright > 2.5, but did not 
find any further increase in the S -statistics. 

B.3. Comparison to Photometric Redshift Surveys 

In order to check the significance of the optimal values em- 
pirically found for mbright and mf a ; nt - i.e., do they provide an 
effective distinction between galaxies at redshift z < 0.5 and 
those at z > 0.5? - and to estimate the geometric factor needed 
to convert gravitational shear into a mass estimate, we com- 
pare our data to two catalogues with kno wn photometric red shift 
distributions, the CFH TLS Deep 1 field (lllbert et al.ll2006h and 
the COSMOS survey (lllbert et al.ll2009l) . In Fig. E31 we com- 
pare the source number counts as a function of magnitude of 
the MMT/Megacam catalogue of the CL0030+2618 field (be- 
fore and after selection of high-quality shape objects, i.e. the 
unflagged SExctractor objects compared to the galaxy shape 
catalogue) to the CFHTLS Dl (MegaCam at CFHT, SDSS fil- 
ter system) and COSMOS photo-z sources. For the latter, the 
Subaru g + r + i + magnitudes similar to the SDSS filters have been 
used. From the CFHTLS we use all unflagged sources classi- 
fied as galaxies, detected in all five bands (u*g'r'i'z') and with a 
photo-z derived from at least three bands whose lcr error mar- 
gin Azp n satisfies Az pn /(1 + z p h) <0.25. Likewise, we use all un- 
flagged sources classified as galaxies having an unflagged photo- 



21 We prefer mf a j nt = 22.5 over the slightly better = 22.4 be- 
cause of the greater robustness of the m faint = 22.5 cases with respect 
to changes in m^ght^ 



z estimate in the COSMOS catalogue that are detected in the 
Subaru g + r + t and CFHT i' passbands. 

Fig ■ [BT4] illustrates how the various cuts in the KSB pipeline 
remove faint galaxies from the catalogue, shifting the maxi- 
mum r' mh of the histogramme from r' mh = 26.0 ± 0.5 to = 
25.0 ± 0.5. We note that the CFHTLS Dl shows a very sim- 
ilar histogramme over most of the relevant magnitude range 
20.5 < r' < 27.0, also peaking at r'^ = 25.0 ± 0.5. The other 
fields of the CFHTLS Deep Survey, D2 to D4, show a behaviour 
similar to Dl and are omitted from Fig. lB.4l for the sake of clar- 
ity. The COSMOS photo-z catalogue, on the other hand, is shal- 
lower, with r\ =24.0+0.5, but showing a number count function 
similar to the one in the CL0030+2618 data at the bright end. 
Therefore, we use CFHTLS as a reference survey, estimating 
the relations between galaxy colours and photometric redshift in 
the CL0030+2618 data from the Dl field and using all fields to 
derive the redshift distribution. 

B.3.1. Photometric Properties 

First, we want to investigate the effect of the photometric cuts 
optimising the aperture mass detection on the redshift distribu- 
tion of the CFHTLS Dl catalogue. 

In Fig. IB.5I we compare the r' - i' and g' - r' colours of 
CFHTLS Dl galaxies with photometric redshift 0.01 <z p h < 0.5 
(upper panel) and z p h > 0.5 (lower panel) to the polygonal re- 
gions found from Fig. IB. 21 containing all bright (/ < 20.0) 
and most of the intermediate (20.0 < r' < 22.5) galaxies in 
the CL0030+2618 field. The bright and intermediate nearby 
(0.01 < z p h < 0.5) galaxies indeed populate a similar region in 
the colour-colour diagramme like their Megacam counterparts, 
albeit being slightly shifted towards bluer g' — r 1 colours. Thus, 
given its simplicity, our background selection is quite efficient 
for the r' < 22.5 foreground galaxies, removing 85 % of them 
from the CFHTLS Dl catalogue. On the other hand, the number 
of bright (f <20.0) background (z p h >0.5) galaxies is negligible. 
Only 28 % of the intermediary CFHTLS Dl background galax- 
ies, redder in r'-i' than the foreground sources but not in g'-r', 
are removed by the selection criteria. 

Concerning the faint (r' >22.5) galaxy population, the first 
observe that, despite the similar source counts (Fig. IB.4I >, the 
colour distributions of faint sources in the CFHTLS Dl and 
CL0030+2618 fields differ qualitatively. Further investigations 
will be needed to relate this observation to a possible cause in 
the data reduction pipeline. This difference in the colour distri- 
butions affects the impact of the background selection: In con- 
trast to the 6.0 % sources removed as foregrounds from the 
CL0030+2618 catalogue, the size of the CFHTLS Dl catalogue 
is reduced by only 0.8 %. (The rates differ little for the D2 to D4 
fields.) 

Second, we note the existence of a significant fraction of 
Z p h < 0.5 galaxies even to very faint magnitudes: we find 15 % 
of the r' > 22.5 sources and 8 % of the r' > 25.0 sources to 
be foregrounds to CL0030+2618, judging from their photo-zs. 
Consequently, our background selection cannot identify these 
sources, leading to a contamina tion of the l e nsing catalogue and 
a di lution of the lensin g signal, lllbert et alJ (120061 their Fig. 16) 
an d lllbert et al] d2009L their Fig. 14) confirm the existence of this 
population of faint galaxies at low z p h. Although there certainly 
is a contribution by catastrophic outliers to which a z p h < 0.5 has 
been assi gned erroneously , the comparison with spectroscopic 
redshifts (Ilbert et al. 2006, their Fig. 12) indicates that most are 
indeed faint nearby galaxies. 
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Fig. B.5. Colour-colour diagrammes of photo-z galaxies in the CFHTLS Dl field. Shown are the r 1 — i' against g' — r 1 colours for 
foreground (0.01 <z<0.5, upper panel) and background (z>0.5, lower panel) galaxies, divided into the three magnitude bins defined 
in Sect. IB.2l and Fig. IB. 21 r' < m^ght = 20.0 (large symbols), mbright <r' < 'Wfaint (medium-sized symbols) and r' > m^m = 22.5 (small 
symbols). Also shown are the polygonals regions giving the locus of bright and intermediate galaxies in the CL0030+2618 field. 
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Table B.2. Best fit parameters zo, A (fixed), and B (fixed) of 
Eq. ( fB~TT ) to the CFHTLS Dl to D4 redshift distributions. 



Field zq A B median(z ph ) (£> ds /D s ) 

Dl 0.87 (1.15) (1.5) 0.91 0.345 

D2 0.76 (1.15) (1.5) 0.79 0.297 

D3 0.80 (1.15) (1.5) 0.83 0.316 

D4 0.90 (1.15) (1.5) 0.95 0.358 



ratios obtained for the four CFHTLS fields (see Table IB.2I > as 
fiducial value and uncertainty margin for the distance factor of 
our Megacam background sample. These values are consistent 
with the results for (Ais/As) computed directly from the his- 
togrammes in Fig. IB. 61 



Hence, applying the background selection to the whole cat- 
alogue, the rate of z p h < 0.5 galaxies only drops from 18.2 % to 
17.6 %. This indicates a similar level of residual contamination 
in the CL0030+2618 background catalogue (Sect. 16.2b . given 
its redshift distribution follows the one in CFHTLS Dl. We ac- 
count for the shear dilution caused by foreground galaxies as a 
source of systematic uncertainty. To this end, we measure 18.0 % 
galaxies at z p h <0.5 in the background-selected CFHTLS Dl cat- 
alogue, once the 2.2 <g'—i' < 3.0 sources, already covered in the 
correction factor for cluster galaxies (Sect. I6.21 > are excised. We 
consider this 18.0 % uncertainty in the systematic error derived 
from shear calibration effects (Sect. |6T6| >. 



B.3.2. Redshift Distribution 

We use the redshift distribution in the CFHTLS Deep Fields to 
estimate (Dd s /D s ), the catalogue average of the ratio of angular 
diameter distances between deflector and source and source and 
observer. In the absence of (spectroscopic or photometric) red- 
shifts of the individual galaxies, this essential quantity has to be 
determined from fields with a known redshift distribution. 

In Fig. IB. 61 we show the binned photometric redshift dis- 
tributions we find for the CFHTLS D 1 to D4 fields after having 
applied the same photometric cuts as to the CL0030+2618 data. 
The apparent spikes seen at certain redshifts in all the four 
fields are artifacts caused by the photo-z determination. Because 
of those, we prefer calculating (D^ s /D s ) using a fit to the 
Znh-distribution. We choose a functional form introduced by 
Van Waerbeke etail l200lb : 
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where zo is a typical redshift of the sources, and A and B shape 
parameters governing the low-redshift regime and the expo- 
nential drop-off at high redshifts. The prefactor including the 
Gamma function renders p z (z v h) a normalised probability dis- 
tribution. We fit the binned redshift distributions in the range 
< Zph < 4, fixing B = 1. 5 for reasons of robustness to the 
default value suggested by IVan Waerbeke et akl (l200lh . Next, 
A = 1.15 is fixed too, to the value preferred for three of the 
four fields. The final results are summarised in Table IB~2l Note 
that -DdsfeUd = 0.5)/D s (z s ) varies substantially over the range 
0.8 < median(zph) < 1 .0 spanned by the median redshifts of the 
fits to Dl to D4 (see Fig. IB.61 >. We now calculate the average 
distance ratio for each field by integrating this function with the 
redshift distribution over all redshifts larger than za = 0.5: 



, ,Ais(z,Zd) 
Pz(z) _ , . dz (B.2) 
D s (z) 



For the mass estimation of CL0030+2618, we use the average 
and standard deviation (Dds/^s) = 0.33 ± 0.03 of the distance 
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Fig. 11. The tangential shear profile of CL0030+2618, averaged 
in bins of 90" width (solid line with dots). Upper panel: The 
best fitting NFW model in the fiducial case (see text; dashed 
line) and the binned cross-component g x of the measured shear 
(dotted line with diamonds). Error bars give the standard devi- 
ation of measured values in the resp. bin. Middle panel: NFW 
models with r2oo from the fiducial fit and concentration set to 
cnfw = 4.0 (dot-dashed line) and fitting to ^200 on ly keeping 
cnfw = 4.0 fixed (triple dot-dashed line), compared to the fidu- 
cial fit (dashed line). Lower panel: The same models as in the 
middle panel, but all truncated at r2oo, with a g t cc 0~ 2 drop-off at 
larger radii. 
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Fig.A.l. Mean stellar anisotropies (ei,2) found in the r'-band 
frames of CL0030+26 18 fulfilling the seeing condition s < 1 .0". 
All images with (|e|) < 0.05 are included in the final coaddition 
(inner circle), while those exceeding (|e|) < 0.06 (outer circle) 
are always rejected. The decision for intermediary objects (see 
below) is based on visual inspection. 
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Fig.A.3. Comparison of the SDSS and Megacam filter sys- 
tems. The plot shows the complete transmission curves for the 
u'g'r'i'z' filters of both systems as a function of wavelength, in- 
cluding the atmospheric transmissivity (as given for the SDSS 
site), the CCD quantum efficiency, and the actual effect of the 
filter, as measured in the laboratory. The solid lines give sensi- 
tivities of Megacam filters for photons incident on the optical 
axis while the dash-dotted lines show the same quantity near 
the corner of the Megacam array. Overplotted as dashed lines 
are the transmission curves defining the SDSS bandpass sys- 
tem. The black, dotted curve shows the Megacam quantum ef- 
ficiency we derive from the instrument specifications, scaled 
by one half to show it conveniently on the plot. Note that we 
need to interpolate its values from only five points in the range 
300 nm</l< 1000 nm and have to extrapolate outside it. 
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Fig. B.l. Upper plot: Colour-magnitude diagramme of KSB 
galaxies with a radial distance 8 < 4' from the centre of 
CL0030+2618. Symbol sizes and shades of grey denote galax- 
ies from the galaxy shape catalog ue in different cluste r-centric 
radial bins. The g'-i' colours of IColeman et alj d!980l CWW) 
template galaxies at z = 0.5 (solid line and large dot at i' = 20) 
and z = 0.25 (dashed line and large dot) are shown for compar- 
ison, as well as four notable bright galaxies detailed in Table [2] 
Lower plot: Colour-colour diagramme with the same objects. 
The polygonal region delineating the locus of bright galaxies (cf. 
Fig. IB.2t is given for comparison. 



Fig. B.4. Source number counts in the CL0030+2618 and exem- 
plary photometric redshift fields. Given are numbers of sources 
as fractions of the total number of objects in the r'-band cata- 
logue for the MMT Megacam CL0030+2618 raw (long-dashed 
curve) and lensing (before background selection; solid curve) 
catalogues as well as for the CFHTLS Dl field (dash-dotted 
curve). The dashed curve denotes the COSMOS r + -band number 
counts. Vertical dotted lines indicate Wbrighi and ;«f a i nt . 
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Fig.B.6. Photometric redshift distributions of the CFHTLS Dl 
to D4 fields after application o f the photometric cu t s defin ed in 
Sect. IB.2I (histogrammes) and IVan Waerbeke et all d2001l) best 
fits to these (solid lines). The function Ddsfc; za - 0.5)/D s (z s ) is 
denoted by dashed lines. 



